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Executive Summary 



The document at hand describes the physics case, the key experiments proposed 
and the technical solution for the project "WASA at COSY" , which concerns the 
transfer of the WASA pellet-target and detector system from CELSIUS (TSL, 
Uppsala, Sweden) to COSY (FZJ, Jiilich, Germany) and its installation at an 
internal target position. 

WASA at COSY will provide unique scientific possibilities for research in hadron 
physics with hadronic probes, since it combines: 

• COSY with its (phase-space cooled, polarized) proton and dcuteron beams 
with energies sufficient to cover the strange quark sector including 
mesons 

• WASA, a close to Air detector for both photons and charged particles 

Since both COSY and WASA are operational devices, the time needed to 
install and commission it is short: it is anticipated that the physics program 
will start at the end of 2006. 

The physics that will be investigated comprises: 

Symmetries and symmetry breaking: 

The fundamental theories of physics are based on symmetry, yet: our world 
is filled with asymmetry (T.D. Lee). 

Finding and studying symmetries and symmetry violations in hadronic 
reactions will help to better understand the strong interaction. Hadron 
beams are particularly well suited for such investigations due to the spin- 
and isospin filter mechanism. In addition, it is even conceivable to search 
for physics beyond the Standard Model. 

Hadron structure and interactions: 

That [intermediate distance] scale is the richest phenomenologically, and 
is certainly the crux region to understand . . . what QCD is really about. 
And at the heart of the subject is the hadron spectrum, in particular the 
spectrum built from light quarks. (■ ■ ■ ) Without question, there is a great 
need . . . for a new round of experiments, especially utilizing hadron beams 
(J.D. Bjorken). 

Finding and further investigating specific hadronic bound systems will — 
together with corresponding progress in theory — provide a more funda- 
mental insight into how nature makes hadrons. 

The proposal concentrates on the study of symmetry breaking as the primary 
objective, in particular in r\ and r/-decays and meson production. In addition, 
the potential of further studies is discussed. 

An experienced and enthusiastic community of more than one hundred scientists 
from all over Europe and abroad has started to work for "WASA at COSY" to 
become a reality and is looking forward to exploit its scientific potential. 
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1 Introduction 



Hadron physics is concerned with one of the most difficult but also most fas- 
cinating problems in contemporary physics. It is commonly accepted that the 
underlying theory, Quantum Chromo Dynamics (QCD), is correct. However, 
very little is known about the solution of this non linear quantum field theory 
in the regime of small and moderate energy scales pertinent to the matter sur- 
rounding us. QCD is characterized by two particular phenomena, confinement 
and chiral symmetry breaking. QCD is formulated in terms of colored quarks 
and gluons. These fundamental particles have never been observed as free states 
but only in composite systems, so called hadrons, i.e. they are confined within 
colorless mesons, baryons and, possibly, exotics. Furthermore, in the sector of 
the light up, down and strange quarks, QCD possesses an approximate chiral 
symmetry, which is not observed in nature but rather broken spontaneously 
with the appearance of almost masslcss Goldstone bosons. These can be iden- 
tified with the lowest octet of pseudoscalar mesons. Such broken symmetries 
also play an important role in phenomenon of superconductivity or the Higgs 
mechanism, which is believed to lead to the creation masses of the quarks and 
leptons. The consequences of the dynamical and also the explicit chiral symme- 
try violation through the current quark masses in QCD can be systematically 
explored with effective field theory. Presumably related to these phenomena 
is the generation of hadron masses in QCD, e.g. the proton and the neutron 
are much heavier than the current quarks they are made of. In the energy 
regime considered here, the interaction between the quarks and gluons is highly 
non-linear and also non-pcrturbative, excluding the use of standard perturba- 
tive methods as applied e.g. in Quantum Electro Dynamics or QCD at high 
energies, where asymptotic freedom is the dominant feature of the theory. A 
deeper understanding of the structure and dynamics of hadrons, in particular 
their excitation spectrum, as derived from QCD would explore one of the last 
terrae incognitae of the Standard Model. 

To understand complex systems such as hadrons, a variety of complementary 
experimental and theoretical studies are mandatory. The WASA detector at 
COSY offers a unique possibility to deepen our understanding of aspects con- 
nected with QCD in the non-perturbative regime through a precise study of 
symmetry breakings and very specific investigations of hadron structure. The 77 
and 77' decays, that vanish in the limit of equal quark masses, allow us to explore 
the explicit isospin symmetry breaking in QCD. 

It is well established that isospin violating quark mass effects are often masked 
by electromagnetic corrections and thus such studies require the detection of 
neutral decay particles. A textbook example for this is the threshold photo- 
production of neutral pions off protons which has so clearly revealed the chiral 
loops inherent to the spontaneous chiral symmetry breaking of QCD. Through 
the same quark mass term of QCD, such quark mass effects also appear in sys- 
tems of a few nucleons, so that dd collisions producing 4 He plus neutral Gold- 
stone bosons at COSY energies offer yet another tool to explore this area - also 
because only now theoretical tools are becoming available that allow a model- 
independent analysis of these processes. Furthermore, precision measuremetns 
of rare 77 and rf decays can be used to get new limits on the breaking of the 
fundamental C, P and T symmetries, or combinations thereof. In addition, 
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reactions involving the rf might also offer insight into the elusive glue through 
the anomaly. 

The spontaneous chiral symmetry breaking is believed to be at the heart of 
the mass generation in QCD. Still, this issue cannot be separated from the 
question: what kind of hadrons are really generated in QCD? In the course of 
the last year, the paradigm of the constituent quark model, i.e. the hadrons 
are either quark-antiquark or three quark states, has been challenged by many 
experiments. However, a definite picture concerning the nature of e.g. the 
light scalar mesons, the ao(980), /o(980), etc., the A(1405) or the pentaquark, 
the O + (1540), has not yet emerged. Evidently, if one docs not understand 
the nature of these states and others, a true comprehension of confinement 
can not be achieved. WASA at COSY is capable to contribute significantly to 
test the models which are offered to explain these states through the precise 
measurements of decay chains or the coupling to other hadrons. It should be 
stressed that all these proposed experiments will be precision measurements, 
because only accurate data will allow to test and constrain the theory. With 
the advent of effective field theory methods, supplemented by coupled channel 
analyses, the theoretical tools of the required precision arc available to analyze 
the data in a model-independent way, reflecting the change of paradigm in 
the foundations in this area of theoretical physics, that in the past has been 
dominated by model building. Moreover, lattice gauge theory offers prospects 
to confront certain observables, in particular spectroscopic issues, with ab initio 
calculations. A very close contact between experiment and theory is foreseen 
and indispensable to achieve progress in the field of hadron physics. 
In this document it is proposed to perform high precision measurements of 
several hadronic reactions to confront the theoretical predictions sketched above. 
In particular, investigations of symmetry breaking in the decay of the 77 and 77' 
mesons and isospin violation in the reaction dd — > air are seen to be very 
promising for this purpose. The combination of the detector system WASA 
and the COSY accelerator facility perfectly matches the requirements needed 
to measure these reactions. The beam energy range, the phase space cooling 
and the availability of (polarized) proton and deuteron beams at COSY are 
essential aspects for high luminosity measurements with low background of these 
processes. The WASA detector provides nearly full solid angle coverage for both 
charged and neutral particles, which is needed to for kinematically complete 
measurements of the multiparticle final states. Furthermore, the usage of frozen 
pellets of liquid hydrogen and deuterium as the target minimizes background 
from secondary reactions but at the same time allows high luminosity conditions. 
Since the dcvelopement and construction of a detector system fulfilling these 
requirements requires a number of years, the transfer of an existing detector 
like WASA to COSY for these measurements will greatly reduce the lead-time 
until physics output can be expected. 

To summarize: The challenge of hadron physics is not only to gain an under- 
standing of QCD in the non-perturbative regime, e.g. to explain the mechanism 
underlying confinement, but also the possibility of unhinging the foundations 
of the Standard Model at low energies. The understanding of the structure of 
matter is intimately linked to the structure and dynamics of hadrons. Thus, 
hadron physics plays a central role in the whole building of physics. This is the 
area of research where WASA at COSY can contribute significantly as will be 
demonstrated in the following sections. 
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2 Physics case 



2.1 Primary objective: 

Symmetries and symmetry breaking 

The degrees of freedom within the Standard Model are leptons and quarks 
for the matter fields and gluons, photons and the weak gauge bosons. How- 
ever, at low and intermediate energies only the leptons and photons play a role 
as dynamical degrees of freedom: and Z° are too heavy and quarks and 
gluons are trapped inside composite objects (hadrons, glueballs, and possibly 
exotics) . The mechanisms behind this process are confinement and spontaneous 
breakdown of chiral symmetry. However, through its symmetries the underlying 
theory is still visible. If a symmetry is spontaneously broken, there is no obvious 
connection between the excitation spectrum and the symmetry. The currents 
are still conserved and, as a consequence, the interactions of the hadrons at low 
energies are largely constrained. The same holds true in the presence of a small 
explicit symmetry breaking: in nature the chiral symmetry is broken by the 
non- vanishing quark masses. Wc will discuss in the following examples, where 
a study of symmetries and symmetry breaking patterns in hadronic systems 
allows to get insights on symmetries and symmetry breaking patterns of QCD. 



2.1.1 Decays of T] and rj mesons 

Isospin symmetry forbidden decays rj(rj') — > 3-7T are sensitive to isospin symme- 
try breaking (ISB) due to the light quark mass difference Am = m^ — m u , and 
provide an approach complementary to ISB estimates from pscudoscalar meson 
masses and production processes (see section 2.1.2 and [1-3]). ISB effects in 
these decay channels can be associated directly with Am, since electromagnetic 
effects are expected to be small [4,5]. Moreover, the ratio of symmetry breaking 
and allowed rj decays to 37r and r\ 2tt, respectively, is directly proportional to 
the square of Am (see Ref. [6] and section 3.1): 

Mixing between isospin eigenstates | 7r° > and | fj > occurs due to the non- 
vanishing matrix element of the quark mass contribution to the QCD Hamilto- 
nian density Ti rn = m u uu + m^dd + m s ss 

< tt° | Km | fj >= A= ((uu - dd) | Um I (uu + dd- ss)) = — ^Am , (1) 

where SU(3) meson wave functions and a pscudoscalar singlet-octet mixing 
angle Q ps f» —19.5° have been employed [7). The mixing angle for the 
physical states ir° and r\ is related to Am by [6,8,9]: 

cos0 TI) sinGv,, \ f tt° \ . ^ y/3Ai 



) = [ n n : ; sin 6^ = ' (2) 

r) J \ -sine,, cos6 m) J \ V J 1 4(m s -m) v ; 

with the average light quark mass m = (m u + m.d)/2. 

Measuring the ratios of isospin symmetry breaking rj — > 3n decays with respect 
to the allowed rj 2ir modes 

_ r (r/ -> AV) _ r(y^A+ 7 r-) 

1 T(r)' -> 777r°7r°) ' 2 r (77' -> tj-k+t:-) ' [) 
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the n-rj mixing angle in this approach can be determined in 77' decays from 

Hi = Pi sin 2 9^, (4) 
with Pi denoting phase space factors [6]. 

The energy dependence of hadronic decays of 77 and 77' to pseudoscalars, i.e. 
77 (77') — > 3 7r and 77' — ► tjtttt, is contained in the Dalitz plot distribution of the 
decay products. Some of the published results obtained for slope parameters of 
the Dalitz plots for 77 — > 3tt [10-17] are not in agreement with each other and 
with theoretical predictions [18,19]. In particular, the predicted slope param- 
eter for the neutral decay is excluded by the most recent published data [17]. 
Preliminary results on both charged and neutral 3-7T decay modes from a high 
statistics measurement at the KLOE facility have recently been reported [20]. 
The energy dependence of hadronic decays 77' — ^ 77 2tt and 77' — ► 37r is discussed 
in the ChPT framework in more detail in [19,21]. The presently available data 
with highest statistics on the 77 2tt decay mode are still inconclusive, whether 
Dalitz plot distributions deviate from phase space. For the 77' — ► 37r decay, no 
experimental data on the Dalitz plot parameters are presently available. 
Hadronic decay modes also constrain the interpretation of the scalar meson 
nonet: Isovector ao(980) exchange is found to dominate 77' — s- 77 27r [19,22], and 
the experimental data are used to fix unknown scalar-pscudoscalar-pseudoscalar 
couplings [22]. In comparison to 77 — > Sir decays, exchange diagrams involving 
/o(980) and ao(980) arc expected to have a stronger influence in 77' decays to the 
3tt channel, since the exchange particles are much closer to the mass shell [23]. 
The axial U(l) anomaly prevents the pseudoscalar singlet 770 from being a Gold- 
stone boson. The dominant contribution to the mass of the singlet state arises 
from the divergence of the singlet axial current, that acquires an additional 
term with the gluonic field strength tensor and, consequently, does not van- 
ish in the chiral limit. Since the physical states 77 and 77' are mixtures of the 
octet and singlet fields, 77s and 770, respectively, the 77' cannot be disentangled 
from the 77 and a proper study of 77' physics requires to account for the 77 as 
well. Precision data on decays of 77 and 77' allow to probe a variety of aspects 
of low-energy hadron dynamics. Experimental results can be compared with 
model— independent calculations of chiral perturbation theory (ChPT) with a 
clear theoretical connection to QCD. Two-photon decays in the 77-77' system 
allows ChPT predictions to be tested, particularly its extension to the U(3) 
framework including the axial U(l) anomaly. 

Radiative decays 77 (77') — ► 2j and 77 (77') — > 7r + 7r~7 test parameters of the QCD 
triangle and box anomalies. Two-photon decays 77 (77') — > 77 provide insight into 
the axial U^(l) anomaly of QCD, since they are determined by the pseudoscalar 
singlet and octet couplings to the divergences of the relevant axial-vector cur- 
rents and the singlet-octet mixing angle O ps (see section 3.7 and [24]). The 
latter can be constrained from radiative decays of vector (V) and pseudoscalar 
(P) mesons in V (P) — ► P (V)7 processes, i.e. from 77' by the p7 and W7 de- 
cay modes. A model framework to conclude about the glue content of the 77' 
is discussed in [25] . The two-photon decays of 77 and 77' arc furthermore well 
suited to confirm the number of colors in the low-energy hadronic sector to 
be N c = 3 [26,27]. Within a few weeks, the available statistics for radiative 
77' decays could be increased by two orders of magnitude at the WASA facility 
(Table 1). 
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Theoretically, the anomalous behavior of the effective action under chiral trans- 
formations is reproduced by the Wess-Zumino-Witten term of the effective La- 
grangian [28,29] which yields the major contribution to the decays 77 (77') — ► 77. 
The 7r + 7r~7 mode is dominated by 77' — ► pj. However, the amplitude of the 
non-resonant 77' — > 7r + 7r~7 decay is determined both by parameters of the AVV 
(axialvector-vector-vector) triangle anomaly, that can be derived from radiative 
decays as discussed above, and by the AAAV box anomaly [28,30-34]. On the 
basis of the presently available statistics, non-resonant contributions have not 
been unequivocally extracted (see Ref. [35,36] and section 3.1). 

Dalitz decays 77 (77' ) — > 77* — > with I = e, [X, where a time-like vir- 

tual photon converts to a lepton pair, probe the transition form factor, i.e. 
the electromagnetic properties of 77 and rf in terms of the spatial distribu- 
tion of meson charge with the l + l~ invariant mass corresponding to the four- 
momentum squared of the virtual photon (theoretical approaches for 77 decays 
are discussed in [33,37]). The available data for 77' — > ~fl + l~ consist of 33 ± 7 
/i + /i~7 events [38,39]. Thus, a detailed analysis of the leptonic mass spectrum 
is presently not possible due to the lack of statistics. Significant progress could 
be made using the WASA facility at COSY (section 3.7). 

The decays 77 (77') — > l + l~l + l~ into lepton pairs address decays via two off-shell 
photons and indicate whether double vector meson dominance is realized in 
nature. The coupling to two virtual photons is substantial for the real part of 
the decay amplitude to a single lepton pair, but also an important issue in kaon 
decays, and for the hadronic contribution to the anomalous magnetic moment 
of the muon [40,41]. The corresponding 7/7*7* form factor has neither been 
measured in the time-like nor in the space-like region (see section 3.1), a recent 
theoretical investigation is presented in [33]. 

The dominant mechanism in the leptonic decay 77 — > e + e~, which is forbidden to 
proceed via a single photon intermediate state, is a fourth order electromagnetic 
process with two virtual photons [39,42]. The decay is additionally suppressed 
by helicity factors m e /m, at the r ye + e~ vertices. A very low branching ratio 
5 • 10 -9 [43,44], predicted by the Standard Model, makes the decay sensitive 
to contributions from non-conventional effects that would significantly increase 
the branching ratio (section 3.7). 

The observation of CP or C violating rare decays might hint at effects from 
new physics beyond the Standard Model: The single conversion decay 77 — > 
ir + Tr~e + e~ allows to search for CP violation in flavor-conserving processes be- 
yond the Cabbibo-Kobayashi-Maskawa (CKM) mechanism [46,47], which are 
not constrained by limits on the neutron electric dipole moment. A recently 
proposed CP violating mechanism from an interference between magnetic and 
electric decay amplitudes can induce a sizable linear photon polarization in 
the 77 — > 7r + 7r~7 decay [48,49]. The photon polarization is accessible in the 
associated conversion decay 77 — > 7r + 7r~e + e~ as an asymmetry in the angular 
distribution between the 7r + 7r~ and e + e~ production planes. Similarly, this 
method was used to measure CP violating effects in the K° system in the de- 
cay Kl — > n + Tr~e + e~ [50], with an asymmetry reported as large as 14%. The 
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Table 1: Counting rate estimates including detector efficiencies for 77 (77') decays 
with WASA at COSY based on a luminosity of 10 32 cm -2 s -1 , with a 20 /_tb 
(300 nb) cross section in pp — > ppr)(rf) at 2.250 GeV/c (3.350 GeV/c) beam 
momentum. 





Decay 


Branching fraction 


Existing 


Counting 




mode 


1 »/ 1 tot 


data 


rate 






[45] 


[events] 


[evts /day] 


V 


e + e~7r + 7r~ 


4. +^ 4 7 u -10- 4 


5 


7000 1 


(semi)-leptonic 




< 6.9 -10- 5 




450 2 




e + e~ 


< 7.7 -10- 5 




1/6 3 




7r e + e~ 


< 4. -10" 5 




1/15- 1/2 4 


rf 


7r + 7T~?7 


44.3 ±1.5 % 


8200 


18000 


hadronic 


7r°7r°77 


20.9 ±1.2 % 


5400 


14500 




3 7T° 


1.56±0.26 TO" 3 


130 


145 




7r + 7r~7r° 


< 5. % 




85 5 


rf 


P°7 


29.5 ±1.0 % 


9550 


44000 6 


radiative 


W) 


3.03 ±0.31 % 


160 


1200 




77 


2.12 ±0.14 % 


2767 


17100 


rf 




1.04 ±0.26 TO" 4 


33 


15 


semi-leptonic 




< 9. -10- 4 




45 7 



corresponding asymmetry for 77 — > tt + tt e + e is expected to be at a level of 
10~ 3 -10- 2 (see section 3.1 and [48,49]). 

A variety of decays, e.g. 77 (77') — > 7r°7r°7, 77 (?/) — ► 37r°7, 77 (7/) — > 37, and 
77' — ► e + e~7r° (77) or 77 — > e + e~7r° probe C-invariance, for which only moderate 
experimental limits exist. For the latter decays to the e + e~7r° mode the domi- 
nant Standard Model — i.e. C conserving — mechanism is an intermediate state 
with two photons 7*7*77° resulting in a branching ratio of 0.2 — 1.3 • 10~ 8 [52-54] 
for the 77 decay. If C invariance is violated, the intermediate 7*77° state with 
one virtual photon can contribute, increasing the branching ratio with respect 
to the Standard Model prediction (section 3.7). 

The move of the WASA facility to COSY offers the possibility to extend the 
programs both on 77 decays at CELSIUS and on 77' production dynamics at 
COSY. Until the end of the present data taking period at KLOE in the middle 
of 2005, an integral luminosity of 2000 pb _1 is anticipated (presently 500 pb~ ), 
corresponding to 8 • 10 7 and 4 • 10 5 fully reconstructed 77 and 77' events, respec- 
tively, from (f> — > 7/ (77') 7 events [56]. The KLOE data will significantly improve 
the presently available data set. However, the total anticipated statistics from 
KLOE can be reached with WASA at COSY in a period of 1 to 10 days, accord- 

1 assuming the theoretical estimate of 3 ■ 10 -4 for the branching ratio [51] 
2 assuming r(j? -» e+ e~ e+ e-)/T to t = 2.52 ■ 10" 5 [51] 

3 Standard Model estimate r(?j — > e+e~)/T tot = 5 ■ 10~ 9 from 77 — > fi+ pr [43,44] 
4 Standard Model estimates r(rj -> 7r°e+e )/T to t = 0.2 - 1.3 • 10" 8 [52-54] 
5 assuming T{q' -> TT+n'TT ) ss T(r]' -» 3tt°), see [19] 
6 including rj' — > 7r + 7r~7 

7 assuming r(t/ -* e+e-j)/T tot 3 • 10~ 4 [55] 



6 



ing to the rate estimates given in table 1, which take into account acceptance 
and reconstruction efficiency of the existing experimental setup of the WASA 
facility. Considering these future developments we concentrate on specific rare 
and not-so-rare decays. 

Together with the decay data, new and important information will be gained 
on the ry'-nucleon interaction. Because of the special features of the 77' outlined 
above, this is a system of interest by itself. In addition, it is well known that final 
state interactions that occur subsequent to short range processes are universal. 
As a consequence the Ty'-nucleon interaction is the same in pp — ► pprj' as in 
B — > ri'pX . Thus, the measurement of the former reaction will provide necessary 
information to interpret the latter one that will eventually be measured at the 
B factories and might shed light on the dynamics of charmless B decays [57]. 

2.1.2 Isospin violation in dd — > oltv 

Within the standard model there are only two sources of isospin violation, 
namely the electro-magnetic interaction and the differences in the masses of 
the lightest quarks [1,2]. Especially in situations where we are able to disentan- 
gle these two sources, the observation of isospin violation in hadronic reactions 
is a direct window to quark mass ratios [2,3] — quark masses themselves are 
not directly observable and additional information is necessary to assign a scale 
to these fundamental parameters of the standard model (see e.g. [58]). 
Already in 1977 Weinberg predicted a huge effect (up to 30% difference in the 
scattering lengths for pir° and nir°) of isospin violation in 7r°7V scattering [1]. 
Also the impact of soft photons was studied systematically [59-62] . Since scat- 
tering experiments with neutral pions are not feasible, in Ref. [63] it was sug- 
gested to use NN induced pion production instead 9 . The authors demonstrated 
that a charge symmetry breaking (CSB) 10 irN seagull term (the Weinberg term), 
required by chiral symmetry, should be even more relevant for Af{,(pn — > dir°) — 
the forward-backward asymmetry in pn — > dir° — than ir — r\ mixing, previously 
believed to completely dominate this CSB observable [65]. 
Interest in CSB in hadronic reactions was revived recently with the successful 
completion of experiments for Afb(pn — > dir°) [67] as well as dd — > cot ( [66], see 
Fig. 1) close to the pion threshold. A large collaboration has formed to perform 
the corresponding calculations within effective field theory. First results are 
presented in [68] . The studies show that the relative importance of the various 
charge symmetry breaking effects is very different for the two reactions and a 
consistent investigation of both should help to further disentangle the leading 
CSB matrix elements. 

However, open questions will remain. It is the main result of Ref. [68] that, 
within a plane wave calculation, the Weinberg term is suppressed due to selec- 
tion rules in the reaction dd — > cot . Thus, the lowest order that gives finite 
contributions is next-to-next-to leading order (NNLO). At this order there are 
several diagrams contributing and only a NNNLO contribution will really prove 
that there is a convergence in the series for this particular reaction. However, 
there is one more important test for the approach that will be used for the 

9 The tt°p scattering length might also be measurable in polarized neutral pion photopro- 
duction at threshold [64]. 

10 Charge symmetry is fulfilled, if the amplitudes are invariant under a 180 degree rotation 
in isospin space. CSB is thus a subclass of isospin breaking effects. 
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Figure 1: Experimental missing- 
mass spectra (dd — > aX) from 
Ref. [66] revealing a clean pion peak. 
The smooth curves show a reproduc- 
tion of the data using a Gaussian 
peak and a continuum. The cor- 
responding cross sections are a = 
(12.7±2.2)pb and a = (15.1±3.1)pb 
for the lower and higher energy, re- 
spectively. The data are consistent 
with pure s- waves contributions. 



analysis of the recent CSB experiments, namely, once the parameters are fixed 
the p-waves in dd — ► an can be predicted parameter free to leading and next- 
to-leading order. The IUCF experiment [66] was consistent with purely s-waves 
contributing. Thus, the same experiment at somewhat higher energies (but still 
well below the A region, e.g. at Q m 60 MeV) is urgently called for. 
The Bose symmetry of the initial wave function strongly limits the number of 
allowed partial waves. As a consequence, for s-wave as well as p-wave production 
only a single partial wave is allowed, namely 3 Pq — » s and 5 Di — > p. Thus, they 
do not interfere in any unpolarized measurement. It would therefore be very 
important to measure the unpolarized total cross section as well as polarized 
differential observables. In addition, the analyzing powers are given by the 
imaginary parts of interference terms and thus would provide an important test 
of the initial state interaction, the ingredient to the calculations that is controlled 
to the smallest extent. In a plane wave treatment the analyzing powers would 
vanish identically. 

In a further step, CSB effects involving the A isobar and purely nucleonic ones 
can be disentangled taking data in the A regime as well (Q ~ 160 MeV). This 
will not only be an additional important cross check of the analysis carried out 
previously, but gives in addition quite direct access to matrix elements difficult 
to get at from other experiments. However, in order to study the different 
partial waves polarized differential observables are as important as for the lower 
energy. 

It is important to stress that CSB studies for few-nucleon systems can be based 
on recent developments in effective field theories [69-73], thus allowing for a 
model-independent analysis of the data. 
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2.2 Additional objectives: 

Exotic and cryptoexotic hadron resonances 

The study of hadron resonances and their decays is one important approach 
to obtain information about quark dynamics in the non-perturbative regime of 
QCD. The experiments performed at the proton-antiproton storage ring LEAR 
(Crystal Barrel) and at the new generation of experiments in electron-positron 
colliders (BaBar, Belle, BES, CLEO) have found new mesons which challenge 
the traditional interpretation of mesons as quark-antiquark states [74,75]. In 
particular, new meson candidates with the quantum numbers J PC = 1 ^ can- 
not be generated by qq states [76,77]. Because of this feature these states are 
called exotics. With the recent evidence for the existence of the so called pen- 
taquark + a first baryon candidate was added to this list. In addition, in the 
hadron spectrum there are several states with quantum numbers consistent with 
the quark model. However, some of their other properties arc not. These states 
are called cryptoexotics. Both mesons, like the well known oq and /o, as well 
as baryons like the A(1405), are discussed in the literature as candidates for 
this class of states. However, it is clear that additional experimental informa- 
tion is needed to verify or falsify these conjectures. Hadronic probes provide an 
experimental tool complementary to electromagnetic excitations because of dif- 
ferent spin- and isospin selectivities. In the following we elaborate why WASA 
is particularly suited to investigate those exotic and cryptoexotic resonances. 
In the meson sector the controversial states that are within the COSY en- 
ergy range are the scalar-isoscalar /o(980) and the scalar-isovector ao(980). 
Although QCD can in principle be treated explicitly in the low Mo-mentum- 
transfer regime using lattice techniques [78], those are not yet in the position 
to make quantitative statements about light scalar states (J p =0 + ). Alterna- 
tively, QCD inspired models, which use effective degrees of freedom, arc to be 
used. The constituent quark model is one of the most successful in this re- 
spect (see e.g. Ref. [79-84]). This approach treats the lightest scalar resonances 
ao//o(980) as conventional qq states. However, they have also been identified 
with KK molecules [85-89] or compact qq-qq states [90-92]. It has even been 
suggested that at masses below 1.0 GeV/c 2 a full nonet of 4-quark states might 
exist [93]. 

As in the reaction dd — ► a7r°, with WASA at COSY we can use the isospin 
selectivity in the dd — > aX transition, to model independently access the isospin- 
violating do-/o mixing: the total cross section for the reaction dd — > a(7r°7y) was 
shown to be directly proportional to the square of the mixing matrix element 
with the influence of other isospin violating effects being suppressed by one order 
of magnitude [94] . There is the additional possibility to extract the phase of the 
mixing matrix element from the forward-backward asymmetry in pn — ► da^ / fo . 
Both measurements will reveal independent new information about the structure 
of the light scalar mesons, in particular about their KK -meson content. 
In the baryon sector, several facilities have obtained evidence for a new baryon 
resonance, the + , which has achieved a three-star status in the Particle Data 
Group [45]. If ultimately confirmed with convincing statistics, this resonance 
will have a large impact on the field. The + has positive strangeness which can 
only be generated by an antistrange quark as a building block of the resonance 
and therefore cannot be accommodated by a three quark structure. The width of 
the + is unexpectedly small: while direct determinations of the width can only 
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give upper limits of the order of 10 MeV/c 2 due to limited energy resolution, 
indirect analysis point to a width of less than 1 MeV/c 2 [95-98]. Jaffe and 
Jain have recently pointed out the limitations such a small width imposes on 
the possibilities to generate the 6 + dynamically in the K N channel or as a 
Castillejo-Dalitz-Dyson pole [99]. 

Theorists have pointed out the possible existence of pentaquark states with nar- 
row widths due to the formation of colored clusters about thirty years ago [100]. 
Recent theoretical work has focussed on the diquark degree of freedom [101]. 
Though the clustering of colored diquarks provides a qualitative explanation of 
the small width of the © + , explicit dynamical calculation require fine-tuning of 
the model [102]. Combining the colored diquark degree of freedom with chiral 
symmetry, however, produces hadrons which are stable and do not couple to 
the kaon-nucleon degree of freedom in the chiral limit [103-105]. 
The production of the + via the reaction pd — > pAQ + allows access to a 
pure isospin 1 = combination AO + in the final state. This allows to obtain 
information complementary to the proton-proton reaction which produces the 
+ starting from a pure isospin 1 = 1 state. 

It is well established that if there is one exotic there should be many and it is 
the corresponding multiplet structure that encodes essential information on the 
underlying substructures. 

The iV*(1710) has been assumed to be a member of the pentaquark antidecu- 
plet and predictions for the branching into various decay channels have been 
published [106] awaiting experimental tests. 

The Roper resonance, 7V(1440), has been studied already with WASA at CEL- 
SIUS. The Roper is not a manifestly exotic state, but its broad width and small 
mass have caused many experimental and theoretical investigations of its struc- 
ture. In the Jaffe- Wilczck diquark scenario, the Roper may be related to the 

e+. 

In addition, the A(1405) has been a candidate for a non-triplet quark structure 
and thus for a crypto exotic for a long time. Recent work within unitarized 
chiral perturbation theory suggests that the A(1405) is made of two overlapping 
resonances with different flavor structure. This calls for detailed investigations 
of the decay channels as a function of the excitation energy. 
Recently, the ANKE collaboration has found evidence for a hyperon Y"(1475) 
(see Fig. 5 in section 2.2.2) which is not predicted by any known quark model. 
The PDG gives a one-star status to a surmised £(1480). 

2.2.1 Mixing of the scalar mesons ao//o(980) 

Both, the (isospin 1=1) ao- and the (/=0) /o-resonances can decay into KK, 
whereas in the non-strange sector the decays are into different final states ac- 
cording to their isospin, aJ^(7T ± ?/)/ = i, aQ-^(7r°7/)/ = i and /o^(7r°7r°)/=o or 
("7r + 7r~)/ = o- Thus, only the non-strange decay channels have defined isospin 
and allow one to (model independently) discriminate between the two mesons. 
It is also only by measuring the non-strange decay channels that isospin vio- 
lating effects can be investigated. Such measurements can be carried out with 
WASA at COSY for tt°- or 77-meson identification, while the strange decay chan- 
nels ao/ fo^KsKs might be measured in parallel. Measurements of the KK 
final state with at least one charged kaon have already been performed at COSY 
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Figure 2: Modulus of the leading term of the mixing amplitude A defined in 
Eq. (5). The two kinks occur at the K+K~ (at 987.35 MeV) and the K°K° 
(995.34 MeV) threshold respectively. 




using magnetic spectrometers. The main results of these experiments and their 
implications on the proposed measurements are described below. 
In 1979 Achasov and collaborators [107] pointed out that between the charged 
and neutral kaon thresholds the leading term to the ao-/o mixing amplitude 
is dominated by the unitary cut of the intermediate two-kaon system. It was 
demonstrated, that the leading piece of the ag-fo mixing amplitude can be 
written as 

A = (/ |T|ao) = i9 f0 KK9 a0 KKVs (pk° - Pk+) + O ( p «° ~ p K + \ , (5) 

where the effective coupling constants are defined through T xK ^ = 9 K kP k ' 
This -y/s dependence of A is depicted in Fig. 2. Here, electromagnetic effects 
between the kaons were neglected as in [107]. 

As demonstrated in a recent compilation [108], the values for g ao xR an d 9f Q KK 
range from 0.224 to 0.834 and from 1.31 to 2.84, respectively — depending on 
the Flatte parameterization of measured 7T7t and irrj spectra. It has also been 
shown in Ref. [108] that the coupling of a physical particle to mesons carries 
information about the nature of that particle, as was shown by Weinberg for 
the deuteron case [109]. Based on the existing data, the authors of Ref. [108] 
conclude that both the do and /o have significant mesonic (KK) components, 
however, more quantitative statements require a better knowledge of g ao /f Q KKi 
which could be obtained from an accurate measurement of A. 
In addition, since A (as defined in eq. 5) is essentially the overlap of the wave- 
functions of ao and /o, it is expected that also that part of the mixing amplitude 
that does not stem from kaons will reveal new insights into the structure of the 
light scalar mesons. 

App — ► dX reaction must lead to (-/=1) production, a,pn — > dX interaction is 
not isospin selective and both the etg and /o may be produced (/=0, 1), whereas 
the dd — > aX reaction — neglecting the small isospin violating contributions 
which are the final goal of the proposed experimental program — is a filter 
for the /o (/=0) resonance, since the initial deuterons and the a particle in 
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the final state have isospin 1=0 ("isospin filter"). Since at COSY it is possible 
to manipulate the initial isospin one can thus selectively produce the ao or /o 
resonances and can identify observables that vanish in the absence of isospin 
violation [3,110]. The idea behind the proposed experiments is the same as 
behind recent measurements of isospin violation in the reactions np—>dir° [67] 
and dd—>air° [66]. However, the interpretation of the signal from the scalar 
mesons is much simpler as compared to the pion case. Since the ao and the /o 
are rather narrow overlapping resonances, the ao-/o mixing in the final state is 
enhanced by more than an order of magnitude compared to isospin violation 
in the production operator (i.e. "direct" isospin violating dd—*aao production) 
and should, e.g., give the dominant contribution to the isospin violating effect 
via the reaction chain dd^afo(I=0)-^aaQ(I=l)^a(ir°r]) [94]. 
The dd—>a(7r°r]) reaction seems to be most promising for the extraction of 
isospin violating effects. Any observation of ir rj production in this particu- 
lar channel is a direct indication of isospin violation, however, the cross section 
da /dm will be given by the product of the mixing amplitude A(m) and the 
dd— >a/o production operator. It is therefore compulsory to determine the latter 
in an independent measurement in order to extract the mixing amplitude. A cor- 
responding experiment aiming at the measurement of the dd — > a(K + K~)i—i—q 
cross section is foreseen for winter 2004/05 at ANKE [111]. These data, together 
with the information on the dd—>a(ir°r]) reaction from WAS A will allow one to 
determine A independently. 

In analogy to the measurement of isospin violation in the reaction np—>dir°, it 
has been predicted that the measurement of angular asymmetries (i.e. forward- 
backward asymmetry in the dao c.m.s.) also allows to extract the ao-/o mixing 
amplitude [112-114]. It was stressed in Ref. [113] that — in contrast to the 
np^dn experiment where the forward-backward asymmetry was found to be 
as small as 0.17% [67] — the reaction pn-^dir°r) is subject to a kinematical 
enhancement. As a consequence, the effect is predicted to be significantly larger 
in the ao/ /o case. The numbers range from some 10% [113] to a few factors [112] 
and, thus, should easily be observable. It has been pointed out in Ref. [114] that 
the analyzing power of the reaction pn~*d'K ri also carries information about the 
ao-/o mixing amplitude. This quantity can be measured at COSY as well. 
An experimental program has already been started at COSY which aims 
at exclusive data on ao//o production close to the KK threshold from 
pp [115-117], pn [118], pd [118,119] and dd [111,120] interactions. The reactions 
pp—>ppK + K~ and pd^ 3 Rc K + K~ have been measured at COSY-11 [115,116] 
and MOMO [119], respectively, at excitation energies up to Q = 56 MeV above 
the KK threshold. However, mainly due to the lack of precise angular distri- 
butions, the contribution of the ao//o to KK production remains unclear for 
these reactions. At ANKE, the reaction pp^dK + K° has been measured ex- 
clusively (by reconstructing the K° from the measured dK + missing mass) at 
beam momenta of p=3.46 and 3.65 GcV/c (Q=46 and 103 MeV). The differ- 
ential spectra for the lower beam momentum are shown in Fig. 3 [121]. The 
background of misidentified events at ANKE is less than 10% which is crucial 
for the partial-wave analysis. This analysis reveals that the K + K° pairs are 
mainly (83%) produced in a relative 5- wave (dashed line in Fig. 3) , presumably 
via the aj channel [121]. Based on model calculations for the pp^dK + K° re- 
action the authors of Ref. [122] draw the conclusion that for Q<100 MeV KK 
pair production proceeds dominantly via the a^-resonance. Limits on the K°d 
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Figure 3: ANKE data for the reaction p(3.46 GeV/c)p -► dK+K° [121]. The 
shaded areas correspond to the systematic uncertainties of the acceptance cor- 
rection. The dashed (dotted) line corresponds to if + .R~ -production in a relative 
S-(P-) wave and the solid line is the sum of both contributions. For definition 
of the vectors p, q and k in the cms of the reaction pp^dK + K° see right hand 
part of the figure. Angular distributions with respect to the beam direction p 
have to be symmetric around 90° since the two protons in the entrance channel 
are indistinguishable. 



scattering length have been obtained in Ref. [123]. For a further discussion of 
scalar mesons in these reactions see e.g. Ref. [124]. 

Data on the reaction p(3.46 GeV/c)p — > dir + X have been obtained at ANKE 
in parallel to the kaon data. In contrast to the latter, where the spectra are 
almost background free, the pp — ► d-K + rj signal is on top of a huge multi-pion 
background [125]. This demonstrates the need for a photon detector supplying 
active 7r°- or r^-mcson identification. With such a detector one may expect data 
of a quality comparable to what is shown in Fig. 3 also for the non-strange decay 
channels. 



2.2.2 Hyperon resonances 

Compared to the spectrum of nucleon resonances, the excitation modes of hy- 
perons (A,£) are much less known. The lowest excitations listed in Ref. [126] are 
the A(1405) J p = \ isospin singlet and the £(1385) J p = | isospin triplet 
states, well accessible in proton-nucleon collisions at COSY. The study of the 
A(1405) state is particularly interesting since this resonance has still not been 
understood in its nature, its spin and isospin assignment is based on indirect 
arguments only. In the quark model the A(1405) is interpreted as a state with 
orbital excitation. However, three well established quark models with entirely 
different residual interaction (gluon exchange [127], meson exchange [128], and 
t'Hooft instanton induced interaction [129]) have difficulties in reproducing the 
mass of the A(1405). All models [127-129] find a degeneracy of the computed 
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Figure 4: Invariant mass spectra for 7r + S and tt E + as observed in the 
p(j,K + ir)T, reaction at SPring-8/LEPS [144]. 



A(1405) with the J p = | A(1520) resonance in contrast to the observation. A 
similar problem is also manifest in very recent lattice QCD calculations where 
an extrapolation to the physical pion mass results in a too high value for the 
A(1405) mass [130]. 

On the other hand, since long time the A(1405) has been interpreted as a KN 
quasi-bound state [131-133], later on substantiated by the consideration of chiral 
symmetry [134, 135]. Recent work in the framework of unitarized chiral effec- 
tive field theories has pointed out that the structure of the A(1405) is extremely 
sensitive to the breaking of the SU(3) symmetry [136-139]. Rcf. [139] finds 
two poles of the scattering matrix close to the nominal A(1405) resonance, one 
of which couples more strongly to 7rE states and the other one mostly to KN 
states. As an experimentally observable consequence, peak structures with dif- 
ferent invariant mass distributions are expected in photon and hadron induced 
reactions that populate the A(1405) resonance [139]. 

The experimental information available on the A(1405) and £(1385) resonances 
is primarily based on studies of K~ proton collisions in bubble chambers [140- 
143]. In this entrance channel the resonances are only observable through the 
decay of higher-lying A* and S* resonances. 

Recently photon- induced hypcron resonance production with a statistics of « 100 
events in the A(1405)/S(1385) region was observed at SPring-8 [144], showing 
different spectral shapes in £ + 7r~ and E~7r + final states. Based on arguments 
given in Ref. [145], this was interpreted as an indication for a meson-baryon 
nature of the A(1405) [144] (see Fig. 4). 

In view of the limited quality of the existing data and the lack of proton induced 
studies, it is important to investigate the nature of the A(1405) resonance in 
proton-proton collisions by measuring its spectral shape with good resolution 
and statistics. This will allow recent theoretical predictions based on the chiral 
unitarity approach to be tested. In particular, a comparison of the E + 7r~ and 
E~7r + spectral distribution will reveal a possible resonant behavior of the the 
1=1 amplitude in this mass region, whereas the £ 7r° distribution provides 
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Figure 5: K + p missing mass spectrum for the reaction p(3.60 GcV/c)p — > 
pK + XT7T+ measured at ANKE in comparison with a simulation (upper panel) . 
The simulated K + p missing mass spectrum decomposes into background from 
misidentified events and non-resonant £~7r + production, as well as the known 
hyperons £(1385), A(1405) and A(1520) decaying into £~7r+ (lower panel). 



the true shape of the 1 = resonance according to Rcf. [145]. A comparison 
with the shape of the spectral distribution measured in photo-induced reactions 
is predicted to reveal the two-pole structure of the resonance [139]. 
An experimental study of the A(1405) resonance also delivers information on 
the £(1385) resonance as a by-product, since due to their similar pole masses 
and widths the production of these two resonances can only be disentangled 
by the complete identification of the final states populated in their decays 
(A(1405) -> £tt (100%); £(1385) -> Att (88±2%), £tt (12±2%)). This involves 
the detection of both charged and neutral pions, and thus a large acceptance 
detector with charged and neutral particle detection capability like the WASA 
detector. In particular, photon detection is required for the observation of the 
decay channel A(1405) — ► £°7r° — > (A7X77). The identification of this final 
state is particularly important since it is not populated in the decay of the 
£°(1385) resonance and thus provides direct access to the spectral shape of the 
A(1405) resonance. The £°(1385) contribution in the charged £7r final states 
(£ + 7r~,£~7r + ) can be deduced from a measurement of the Air final state which 
is not populated in the A(1405) decay. 

Recent data from ANKE on the reaction pp -» pK+Y* -> pK + £(1190)-tt+ 
measured at a beam momentum of 3.60 GeV/c are consistent with a hyperon 
resonance Y* with a mass of 1475 MeV/c 2 and a width of 45 McV/c 2 (see [146] 
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and Fig. 5). This state is also seen with comparable strength in the charge- 
conjugate decay channel Y* — ► £(1190) + 7r~ and its production cross section 
is of the same order as for the A(1405). Within the above mentioned quark 
models [127-129] there is no room for an additional A or £ hyperon below 
«1600 MeV/c 2 . However, a £(1480) with a iqlq structure has been predicted 
in Rcf. [100] — as the I — 1 partner of the A(1405) in a common multiplct. 
Further studies are needed in order to confirm the conjecture that the observed 
enhancement in the K + p missing mass spectrum originates from a new hyperon 
resonance. If this interpretation will be supported after further analysis or new 
measurements, the next step is to determine the isospin of the new hyperon 
state. This is not possible from the ANKE data but requires a comparison of 
the possible decay channels , £ 7r°, and Att°. 

As an additional note, it should be mentioned that the different nature of the 
A(1405) and £(1385) states is also expected to manifest itself in different in- 
medium properties of these hyperon resonances. In particular, if the meson- 
baryon picture of the A(1405) is correct, its medium properties arc intimately 
related to the behavior of antikaons (K~ ,K°) in nuclear matter, which is a still 
unsolved problem, of importance in both hadron physics and astrophysics. A 
comparative study of A(1405) and £(1385) resonance production on nuclear 
targets allows to investigate the propagation of these hyperon resonances inside 
the nuclear medium. A better knowledge of the in-medium properties of the 
A(1405) (absorption cross section, modification of the spectral shape) will also 
help to understand the dynamics of antikaons in nuclear matter. 



2.2.3 Pentaquarks 

Early in 2003, the LEPS Collaboration at the SPring-8 facility in Japan ob- 
served a sharp resonance, 6+, at (1.54 ± 0.01) GcV/c 2 with a width smaller 
than 25 MeV/c 2 and a statistical significance of 4.6<r in the reaction 771 — ► 
K + K~n [147]. This resonance decays into K + n, hence carries strangeness 
S = +1. Later, many other groups have claimed the observation of this 
state [148-158]. In particular, evidence for the + has also been seen at the 
COSY-TOF experiment in the reaction pp — > J^ + pK s [159] in the pK s invariant 
mass spectrum. While in some experiments detecting a pK s in the final state the 
strangeness of the resonance in this system is undefined, in this experiment it 
is uniquely determined that the observed pK s system originates from a S = +1 
state since the full final state with an associated £ + is identified. 

Since all known baryons with B = +1 carry negative or zero strangeness, such 
a resonance must have a minimum quark content uudds, a structure which 
is clearly beyond the conventional quark model. This experimental discovery 
triggered a lot of theoretical activity and up to now over 250 papers appeared 
trying to interpret this exotic state [160]. 

There is preliminary evidence that the + is an iso-scalar because no enhance- 
ment was observed in the pK + invariant mass distribution [150,152,153,158]. 
All the other quantum numbers including its angular momentum and parity 
remain undetermined. Most of theoretical work postulated its angular momen- 
tum to be one half, because of its low mass, but the possibility of J = 3/2 still 
cannot be excluded completely. 
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It is important to point out that many other experimental groups reported 
negative results [161-163]. A long list of experiments yielding negative results 
(still unpublished at the time of the writing) can be found in Ref. [164]. 
With both positive and null evidence for the + from a variety of experiments, it 
is difficult to conclude whether the + exists or not. Furthermore, the theoreti- 
cal difficulties to explain a possible narrow width, perhaps as small as 1 MeV/c 2 , 
suggest that if the + exists, it is very unusual indeed. Guidance from lattice 
gauge theory is not helpful, as some calculations show evidence for a negative 
parity resonance, while another indicates positive parity. Furthermore, one does 
not see a resonance in either parity. So the confirmation of the + is of high- 
est importance for which the COSY-TOF collaboration will start an extensive 
experimental investigation in October 2004. 

One should note that the + is a member of a multiplet. This calls for the 
existence for additional unusual states. WASA at COSY with its very good 
mass resolution and large acceptance can make a decisive contribution to this 
field by providing high statistics data on the differential cross sections of exotic 
baryon production in elementary reactions. While there are several elementary 
reactions with the + in the final state within the COSY energy range, the 
pd — > pAQ + reaction is of special interest. Currently this reaction is investi- 
gated by the CELSIUS/WASA collaboration close to threshold. This reaction 
provides access to the pure isospin zero 6 + production channel via the subsys- 
tem pn — > + A and thus allows to obtain further and independent insight into 
the physics of the + pentaquark as compared to elementary transitions like 
pp — > + S + (which is pure isospin one). For example, one can obtain additional 
and complementary information on the parity of the + [165] but explore also 
the + production mechanism, which is so far completely unknown. 
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3 Key experiments 



In this chapter several key experiments are presented to address the issues dis- 
cussed in the previous chapter. After the commissioning phase of WAS A at 
COSY, in which the performance of the detector will be demonstrated, the 
intention is to carry out the experiments discussed below. The proposed imme- 
diate program consists of experiments that are feasible shortly after the initial 
commissioning of WASA at COSY, and will produce physics results soon. The 
medium term plans require either extended data taking or improvements of the 
experimental conditions depending on the experience gained in the initial exper- 
iments. The order of the experiments listed below presents a reasonable balance 
between the physics interest and timeliness of the experimental feasibility. 

3.1 Not-so-rare rj' decays 

Isospin symmetry breaking in 77' decays. The experiment aims at a pre- 
cise determination of the ratios of isospin symmetry breaking decays 77' — ^ 3 7r 
with respect to the isospin allowed 77 2w modes as defined in relation (3) of sec- 
tion 2.1.1 to extract the mixing angle O^,,, which is proportional to the light 
quark mass difference Am (eq. 2). While only an upper limit is known for 77' — > 
7r°7r + 7r~ [45], in the neutral decay modes in eq. (3) the ratio IZi has been mea- 
sured with a statistics of « 130 events, leading to sin 9,^ = 0.023 ± 0.002 [166]. 
However, there is considerable uncertainty in the literature concerning sinO^,,: 
a value of sin @ nr] = 0.010 was extracted using PCAC relations for meson masses 
and removing electromagnetic contributions [6] , in agreement with the estimate 
in [167]. Further results given in the literature cover a range between 0.013 and 
0.014 [168-170], except for sin 6^ = 0.034 ± 0.013 reported in [171]. 
With the WASA facility at COSY, in the tagging reaction pp — > ppr] 1 a statistical 
accuracy of 1 % can be achieved for the lowest sin O w)) estimate within a running 
time of 12 weeks, improving the precision of the existing experimental value by 
an order of magnitude (see Table 1). Background influencing the precision of 
the extracted value is primarily expected from non-resonant 3 7r production. 
Although the background will be distributed smoothly in the vicinity of an 77' 
signal in the two proton missing mass, the final precision of the result will depend 
on the four-momentum resolution achieved for both protons in the forward 
detector and the decay system. Modifications of the WASA detector setup to 
ensure a required performance at higher energies are discussed in section 4.3. 
With selective central detector trigger conditions on either six neutral particles 
from the decays 77(71-°) — > 27 or two neutral and two charged particles, both ratios 
IZi can be determined simultaneously Since the Dalitz plots of the hadronic 
decay channels arc measured completely, the energy dependence of the decays 
is inherent to the data, and can be used both to test predictions from ChPT, 
and as a constraint on the the scalar meson sector as discussed in section 2.1.1. 
The present performance of the WASA detector for the detection of multipionic 
final states is demonstrated in Fig. 6 for a data sample obtained at CELSIUS at 
a kinetic beam energy of 1360 MeV, i.e. at an excess energy of 40MeV with 
respect to the 77 production threshold in proton-proton scattering [172]. 
The experimental pp missing mass spectrum is well reproduced by contributions 
from pp — > ppr\ — > pp3ir° and prompt 37r° production (fig. 6a). The complete 
kinematical information is contained in the Dalitz plot for the 77 — > 3 ir° decay, 
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Figure 6: a) Measured pp missing mass 
(crosses) and Monte Carlo event mixture 
(shaded area) of prompt 3 7r° (dashed line) 
and resonant r\ — * 3 7r° production (solid 
line) . b) Symmetrized experimental Dalitz 
plot for the decay rj — > 3 7r°. c) Efficiency 
corrected, normalized radial density distri- 
bution of the Dalitz plot. Crosses denote 
data, the solid line corresponds to the fit of 
the slope parameter a (figures from [172]). 
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that is shown in its symmetrized form (< T T > — T^ vs. T^ — T^i/V^ ) in 
Fig. 6b. The slope parameter a (fig. 6c), which reflects the strong and energy 
dependent tttt interaction, is extracted from a linear fit of the normalized radial 
density distribution [7]. The preliminary value a = —0.03 ± 0.025 is limited by 
the statistics of 11700 events, but agrees in sign with the result a = —0.031(4) 
of the Crystal Ball analysis based on 10 6 events [17]. 



Search for evidence of the box anomaly of QCD in r/' — ► 7r + 7r 7 decays. 
The experimental signature for the detection of the 7r + 7r~7 decay mode of the 
77' consists of two charged and one neutral particle in the central detector, with 
the rj being tagged from the missing mass with respect to the two protons iden- 
tified in the forward detector. This signature can be used for a selective trigger 
condition. Moreover, with the requirement of one 7 in the central detector, the 
experiment trigger for a measurement of the hadronic decays rf — * ?7(7r )7r + 7r~ 
discussed above is a mere subset, i.e. both experiments could run simultaneously. 

The experiment aims at measuring the 7r + 7r~ invariant mass distribution in 
the decay 77' — > 7r + 7r~7 with high accuracy. One (Two) order(s) of magnitude 
higher statistics compared to the presently available data can be achieved within 
a running time of two days (three weeks). Following the experimental technique 
used in [36] the number of 7/ — * 7r + 7r~7 can be determined in each bin of the 
7T + 7r - invariant mass separately. Consequently, the 7r + 7r~ spectrum will be 
background free (Fig. 7). 
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Figure 7: Distribution of the n + ir~ invariant mass for 77' — > 7r + 7r~7 decays 
from [36] . Crosses denote experimental data, the dashed line is the result of the 
fit used to extract box anomaly parameters with a statistical significance of 4 a. 



The quantity E v / associated with the matrix element of a non-resonant contri- 
bution as a consequence of the box anomaly is related at the chiral point to the 
parameters of the pseudoscalar nonet by one of the Chanowitz relations [30,31]: 



E v , (0) 



1 



4tt 2 v / 3./^ 



sin 6 



ps 



^cosOp, 
fo 



(6) 



O pij denotes the pseudoscalar octet-singlet mixing angle, /„. is the pion leptonic 
decay constant, and, in analogy fo and fa are defined as couplings of the pseu- 
doscalar singlet and octet states 770 and r)$ to the divergences of the singlet and 
octet axial-vector currents, respectively. E v > can be derived from a fit to the 
spectrum of the it + it~ mass (Fig. 7 and [36,173,174]) in the ir + ir~j decay 
mode of the 77' using (for a more sophisticated analysis of this decay, see [34], 
which should be used in the final analysis) 
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where k 7 and q,,. are the four-momenta of the outgoing 7 and 7T , g I? / p7 denotes 
the coupling constant for the p contribution to the 7r + 7r _ 7 decay mode, and 
Dp (m W7r ) and G p (m W7r ) are the p propagator and its coupling to the 7r + 7r~ 
channel. 

On the one hand, the extraction of the box anomaly contribution is strongly 
dependent on the model employed in the description of the p. On the other 
hand, from the two-photon widths of 77 and 77' as described by the Chanowitz 
equations [30, 31] using additional information from radiative J/ty — > 77(77^7 
decays (AFN relation [175,176]), the parameters of the pseudoscalar nonet are 
completely determined, and can be used to predict the box anomaly constant 
E v > . Consequently, the p shape and the box anomaly parameters are closely 
related to each other — the approach might allow for an alternative approach 
to extract the p line shape in a clean way. 

So far, experimentally no consistent picture has been obtained: The box 
anomaly contribution extracted from the Crystal Barrel data [36] is rather 
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Figure 8: ir + ir~ effective mass distribution in the decay rj' — > ir + ir~j from 
L3 [35]. Points denote experimental data, the solid line is the best fit result 
with m p = 766 MeV and T p = 150 MeV. The estimated background (shaded 
area) is mainly attributed to decays of the 02(1320). 



disfavoured (confidence level 3%) by the L3 result (Fig. 8 and [35]), which 
is consistent with a pure p line shape (C.L. 37%). On the other hand, the 
Crystal Barrel data would require a mass m p = 790 MeV when neglecting the 
anomaly related non-resonant contribution, i.e. 20 MeV above the presently ac- 
cepted value [45]. In this respect, these two results with the up to now highest 
statistics arc inconsistent. 



3.2 a+(980) production in pp^-dir + rj 

a^-production has first been investigated at COSY by measuring the reaction 
pp — > dK + K° at a beam momentum of p = 3.46 GeV/c (Q = 46 MeV) with 
ANKE [121]. The measured total cross section for the pp — > dK + K° reaction 
is (38 ± 2 ± 14) nb and, thus, wlOOO events could be collected within five days 
of beam time using a cluster-jet target (L — 2.7 • 10 31 cm _2 s~ 1 during these 
measurements). According to the partial- wave analysis the K + K° pairs are 
mainly (about 83%) produced in a relative S- wave, i.e. via the ajj" channel [121]. 
This observation is in good agreement with a model prediction [122] for the total 
a^-production cross section in pp collisions. This model has also been used to 
estimate the cross-section ratio for resonant (via the aj) [112, 177] and non- 
resonant [112,178] ir + r) production: 

pres/nres | o{PP ~> da+) 

R d | P =3.46GcV/c ~ a{pp ^ d7T+v) ~ °- 3 ' • ■ °- 5 - (8) 

This prediction for jj™ s / mcs j s m K ne w jth data from ANKE for the reaction 
pp — > dir + X [179]. Since most of the non-resonant ir + r) pairs are expected 
at lower invariant masses, we anticipate that the resonant signal can well be 
identified in case of the pp — > dn + ri reaction. 
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Figure 9: Invariant mass (ir + r)) for production in pp^da^ — >d7r + r\ and for 
the two background processes, pp — > dn + rj and — > pnn + ri (downscaled by 
factor 10). 



The model can also be used to estimate i?p^ nrcs for the pp — > pmr + r] reaction 
(see Refs. [180,181] for Oq- and Ref. [178] for non-resonant n + r] production): 

prcs/nres I _ a (PP ~^ P na ) _ n m r n n o f Q \ 
H ' p=3.46GcV/c — — T T ~ 0.015 . . .0.1)3. (9) 

' IK ' <J(PP — * pmr+ij) 

Thus R T pn / nrcs is expected to be about one order of magnitude smaller than the 
corresponding ratio ^ es / nres w Jth deuteron formation in the final state. There- 
fore, an experimental study of the in pp reactions with WASA requires the 
identification of the pp — > dn + ri reaction, and the pp — > pnn + r] reaction will be 
the most "dangerous" source of background. Consequently, the following simu- 
lation calculations focus on the problem of deuteron-vs. -proton discrimination. 
The reaction pp^da^ — >d-7r + r\ with subsequent decay ry— >77 can be measured 
with WASA by detecting deuterons in the forward detector (FD) in coincidence 
with the 7r + and photons in the central detector (CD). The reaction is identified 
by reconstructing the masses m(r\) = 771(77) an d m(d) = m.m. (7r + 77). In order 
to investigate the acceptance and background suppression the simulations were 
performed for aj(980) production with a deuteron in the final state and both 
background reactions using the cross section ratios a(aQ):a(dn + r]):a(pmT + 7]) = 
1.1 /jb:3.5 /J,b:96 /ib [112,177,178]. Fig. 9 shows the initial [Tr + r]) invariant mass 
distributions for these reactions: a Flatte distribution for the a^(980) [177,181], 
a distribution according to the model in Refs. [178,181,182] for the non-resonant 
dit + ri, and a phase-space distribution for the non-resonant pnir + rj production. 
The WASA acceptance for forward going particles is 9 w 3° . . . 18°. This accep- 
tance covers about 93% of deuterons from pp^da^ , which are distributed in the 
range = 0°... 17°. The acceptance for pions and photons is 9 w 20° . . . 169° 
(for the simulations angles 9 « 20° . . . 143° — SEC and SEF, the central and 
forward part of the calorimeter — were used). The result of the acceptance 
estimate is shown in the Table 2. For the background reactions the acceptance 
is smaller due to the wider angular distributions. We have only considered the 
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Figure 10: Initial (dotted line) and reconstructed (solid) invariant mass (t: + ^j) 
(left) and missing mass (-7r + 77) (right) for the reaction pp^>dd^ '— >d7r + r\. 



decay channel 77— >77 since for the other two main decay channels 77 — ^7r + 7r tt' 
and 77— >37T° the acceptance is two times smaller. 



Reaction 


Acceptance 




(77^77 decay only) 


pp — >d&Q — >dir^ 77 


0.42 


pp — > dir + ri 


0.34 


pp — ► priTr + i] 


0.32 



Table 2: Simulated acceptance of WASA for a^(980) production and the 
two background processes. Deuterons/protons are detected in the FD (0 f=a 
3° . . . 18°), pions and photons in the CD (0 « 20° . . . 143°). 



In the simulations the ir + detected in the CD have angular resolution better 
than a(0) = 0.4° (for angles larger than 30°) and momentum resolution between 
<?(p)/P = 3 - ■ - 6 %- For photons in the CD the resolution is w 5°(FWHM), due to 
the crystal sizes in the SEC. The energy resolution varies from 10 MeV (FWHM) 
for large angles and low initial energies up to 40 MeV (FWHM) for angles less 
than 30° and initial energies close to 1 GeV. Such resolutions yield 77 invariant 
mass distributions with a FWHM of 43 MeV/c 2 in the 77 mass region, which is 
in agreement with WASA at CELSIUS data [183]. 

Fig. 10 shows the initial and reconstructed invariant mass (7r + 77) and missing 
mass (-7T+77) for the reaction pp—fdaQ—td'K + rj—fd'K + ^"f. Pions and photons 
detected in coincidence in CD were used for the reconstruction, which provide 
a (7r + 77) missing mass resolution of a = 40 MeV/c 2 . 

The expected background from the reaction pp — ► pnix^i] is two order higher 
than the signal. Due to the large momenta of the protons and deuterons 
they cannot be stopped in the FD and their initial kinetic energy cannot be 
reconstructed. Moreover, their energy losses are close to minimum ionizing 
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Figure 11: Reconstructed missing mass (tt + 77) for the reactions 
pp—>da,Q —k!t{ + ?y and pp — ► pmr + rj (left). The cut is indicated by the arrow. 
Energy losses in the FD vs. kinetic energy of the deuterons/protons (calculated 
from the detected ir + and the two photons assuming that the forward particle 
is a dcuteron). The cut is indicated by the two lines. 



and the standard WASA at CELSIUS AE/E method cannot be used for p/d 
discrimination. 

In order to suppress the proton background a set of criteria has been applied. 
The first cut was applied on the reconstructed (7r + 77) missing mass, see Fig. 11. 
The cut at 1.95 GeV suppresses protons by a factor as 1.6. The two-dimensional 
distribution energy loss vs. kinetic energy of the forward particle shows a corre- 
lation for a q events, see Fig. 11. Applying a gate with ±20 MeV around these 
events, background protons can be suppressed by a factor of as 3.9. 
Another strong criterion for proton suppression is the difference between the 
measured azimuthal and polar angles of the forward particles and the expected 
deuteron angles calculated from the tt + and two photons. For the events 
these angles coincide within the resolutions of the detector and the reconstruc- 
tion procedure. For protons from the pmr + r] events the correlation is much 
weaker. This is clearly seen in two dimensional plots of measured vs. calculated 
angles. If one assumes that all forward particles are deuterons, then the real 
deuterons are seen as lines whereas the protons are smeared out (Fig. 12 (a,b)). 
Cuts were applied to the difference between measured and calculated azimuthal 
and polar angles of forward particles (Fig. 12 (c,d)). An azimuthal angle cut of 
±20° and a polar angle cut of ±2° suppresses protons by a factor as 21. 
Taking into account all mentioned cuts and the difference in acceptances the 
proton suppression factor is 1.3x1. 6x3. 9x21«170. We therefore expect that the 
non-resonant pmr + r] events can sufficiently be suppressed without modifications 
of the existing FD for the higher particle momenta at COSY. 
Assuming a luminosity of L = 10 31 cm _2 s _1 , an overall efficiency due to detector 
acceptance, reconstruction algorithms and d/p cuts of as 0.07, the effect of dead- 
time and detector efficiencies of as 0.5 and using the cross section estimate from 
Refs. [112,177], a(pp->da+^dn + ri)xBR{r) -> 2 7 ) = l.l^bx0.393, the final 
count rate is as 0.15 s^ 1 as 90000 week -1 . 
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Figure 12: (a,b) Measured vs. calculated forward azimuthal and polar angles 
assuming that the forward particle is a deuteron. (c,d) Difference between 
measured and calculated angles of the forward particles. The background pnir + r] 
events (dotted line) is downscaled by a factor 100. The cuts are indicated by 
the arrows. 



3.3 Pentaquarks 

At present, several experimental groups, including TOF at COSY, are perform- 
ing high statistics experiments to confirm the existence of the 9 + . 

Provided the existence of the + is confirmed, we will be interested in basic 
properties of this hadron and details of the production mechanism. 

In the following we concentrate on the pd — > pAO + reaction as one example 
to illustrate the feasibility of high statistics experiments involving both neutral 
and charged decay products. 

Out of the two decay modes of the + , the K°p channel gives the best exper- 
imental conditions for being measured in the WASA detector, so we shall here 
concentrate on this decay branch. 
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For the A there are two important decay modes, pn~ and nir° with branching 
ratios of 63.9% and 35.8%, respectively. The choice of using the charged decay 
mode to identify the A is advantageous due to the following reasons: 

• the energy resolution of WASA for charged particles is generally better 
than for photons, 

• the combination of the forward tracker and the MDC provides a clean 
identification of the A by reconstruction of its decay vertex, 

• both decay products are measured, 

• higher branching ratios as compared to neutral decay modes. 

Regarding the kaon decay, half of the events are lost due to K\ having cr = 
15.5 m and thus escaping detection. For the Kg two modes are important, 
ir+ir- {BR = 68.6%) and 7r°7r° {BR = 31.4%). Selecting the charged mode 
has the advantage of higher rate. However, it leads to a final state with six 
charged particles. It might be difficult to distinguish between the ir~ produced 
in the A decay and the Kg decay, which could impair the reconstruction of the 
A. Therefore, the neutral decay mode will be used. 

Taking all branching ratios together gives the overall branching ratio, or effi- 
ciency of detecting 6+ using this particular final state, of 0.5x0.639x0.314 = 
0.1. This final state with four charged particles, 7r~ and three protons, and 
the four photons from tt°tt provides a well defined trigger as well as favorable 
analysis conditions. 

A rough estimate of the WASA acceptance for the pd — > p+A+Q + — > p+pTT~ + 
P7777 reaction has been obtained using phase space Monte Carlo simulation. 
The first things needed are the energies and angles of all charged final state 
particles. They are shown in Fig. 13 below. The proton produced directly is 
denoted p while those produced in the decays of the A and 6 + are denoted pa 
and pq, respectively. It is seen that the maximum energies for p and p\ are 
below 300 MeV, which is the highest energy deposited by protons in the FD. 
For such protons one may expect a very good energy resolution. 
Assuming that reconstruction of the A requires the p\ proton to be detected 
in the FD (2.5° < 9 pA < 18°) and the n~ either in the FD or MDC (22° < 
9 V < 158°) one finds that this requirement is fulfilled in m 85% of cases. This is 
roughly the efficiency of the A reconstruction, e A . It should be noted that the 
resolution of the MDC is limited at the forward and backward angles due to a 
reduced number of wire planes being crossed by the particles. All 17 layers are 
crossed only for angles between 44° and 134°. For the ir~ going into the FD or 
into this limited angular range of the MDC, the £a is reduced to 55%. 
The requirement that the proton pq from + enters the FD or the plastic barrel 
is fulfilled in e p e = 87% of cases. 

Finally, the requirement of the four photons (from the Kg decay into 7r°7r°) to 
fall into the angular range of the electromagnetic calorimeter is fulfilled with an 
efficiency £4 7 = 62%. This gives the overall acceptance of w 45% for eA = 85% 
and 29% for e A = 55%. 

Remembering that the resolution of the WASA detector is far better for charged 
particles than for photons one would like the 9 + mass to be calculated as the 
missing mass to the reconstructed p, p\ and n~ (for 7r _ only angles will be used, 
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Figure 13: Energy and angular ranges of the four final state charged particles 
in the reaction pd — > p + A + + — * p + pjt~ +P7777 at 1360 MeV proton beam 
energy. 



see below) . In such a case one has to take into account losses of protons due to 
nuclear interactions with the detector material, which is expected to reduce the 
acceptance by roughly a factor of two. 

As already mentioned, the + mass will be found as the missing mass calculated 
for the p and p\ reconstructed in the FD and ir~ cither in the FD or MDC. 
Measurement of the energy and angle of the p proton is straightforward since 
the reaction vertex is at the known target position. This is not the case for the 
p\ because of the A's cr=7.89 cm. In order to reconstruct the A momentum one 
has to find the A decay vertex using either the forward tracker, when both pa 
and 7r~ enter the FD, or the forward tracker and the MDC. Reconstruction of 
the vertex requires finding two tracks which, together with the target point, are 
in the same plane. When two tracks fulfilling the above criteria are found, one 
can calculate the absolute value of the tt~ momentum, k, from the condition: 
{yjml + k 2 + E) 2 - (ku + p) 2 = m\, 

where E, p represent the proton total energy and momentum, u is the unit 
vector representing direction of the presumed 7r _ and ma is the A mass. One 
should verify that the reconstructed A-momcntum points to the target area. 
A further check is provided by the AE-E method when the pion enters the FD 
(see Fig. 14) and extracting the momentum from the measurement of the track 
curvature in the magnetic field, when the pion enters the MDC. 
At present the position resolution of the forward tracker amounts to ks 2 cm 
in the x-y plane (perpendicular to the beam) [184] so one can estimate the 
resolution in angle to be of the order of arctan(2/140) « 1°, where 140 cm is 
the distance between the target position and the FD tracker. For the proton 
p, the resolution in angle is given by the size of the beam-target overlap and 



27 



JLJ U ill ' JU ' iJ L.UL L.Jl \LJ l.t. 

Figure 14: AE-E plots for three consecutive layers of the range hodoscope in 
the FD as obtained in the M-C simulation of the reaction pd — > p + A + + — ► 
p + Ptt~ +P7777 at 1360 MeV proton beam energy. 
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Figure 15: Phase space Monte Carlo simulation of the reaction p+d — > p+A+Q + 
at 1360 MeV beam energy. The 9 + peaks are seen on a wide background from 
non-resonant p + d^p + A + K a + p reaction, which has been assumed to be 
10 stronger than the + production. For both peaks the energy resolution for 
protons has been assumed to be equal to 3% and the 7r~ angular resolution to 
amount to 3°. The narrower peak, FWHM=4.6 McV/c 2 , corresponds to the 
proton angular resolution of 1° and for the wider one the resolution of 3° has 
been assumed. The width obtained in this case amounts to 9.6 MeV/c 2 . 



accuracy of the FD tracker. For the ir , which in most cases enters the MDC, 
one can expect the resolution in 9 to be of the order of arctan(l/22) = 2.6°. 
A rough estimate of the resolution in Mq+ can be obtained by smearing the 
energies of the p, p\ and angles of the p, p\ and tt~ with the experimental 
resolutions. Taking for the energy resolution 3%, and 1.0°, and 3.0° for the 
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angular resolution for protons and tt~ respectively, leads to the results shown 
in the Fig. 15 (solid line). The width of the peak at half maximum amounts 
to 4.5 MeV/c 2 , which is nearly five times better than the resolution quoted 
in Ref. [147]. This is a valuable feature which could improve our knowledge 
of the true width of the + resonance expected to be much smaller than the 
20 MeV/c 2 measured in the experiments so far. The wide background seen in 
the figure corresponds to the phase space simulation of the non- resonant reaction 
p + d -> p + A + K Q + p. 

For the + cross section of 30 nb, the beam-target luminosity of 10 31 /cm _2 s _1 
and the acceptance of 1.5% (with the overall branching ratio of 10% included) 
one would expect the 6+ event rate of 30 x 10~ 33 x 10 31 x 0.0 1 5s _1 = 0.0045s" 1 
at a beam energy of 1360 McV. Simulation shows that the detector acceptance 
depends rather weakly on the beam energy. Thus increasing the beam energy 
to 1450 MeV would increase the event rate by a factor of ~7 due to the change 
of the cross section [185]. With a duty factor of 50% obtaining 10 4 events would 
require less than 200 hours. 

3.4 Isospin violation in dd — > air 

Due to the small cross section (a few pb) and the high background from the reac- 
tion dd — > aX the detection and unambiguous identification of the air channel 
is rather difficult. Hence, first attempts to measure this reaction have produced 
only upper limits (see Ref. [186]). The only positive report at a dcutcron energy 
of 1.1 GeV [187] has been questioned, because the background from the double 
radiative capture process dd — > 077 could have been misinterpreted [188]. Con- 
sequently, in the recent experiment performed at IUCF [66] the a and the two 
7's have been measured in coincidence to provide a clean signal of dd — > an 
channel. 

At COSY, studies of this reaction were initially proposed for BIG KARL [189]. 
However, a sufficient background suppression without photon detection would 
have been hard to achieve. The idea was resubmitted within a proposal for an 
electromagnetic calorimeter at ANKE [190]. This combination — like WAS A 
- would have provided a clean identification of the forward-going a particle 
and the ir° decaying into two photons. 

Compared with the latter proposal, the WASA detector provides a similar ac- 
ceptance for photons (O = 20° . . . 169°) but a significantly wider angular accep- 
tance for a particles (O ss 3° . . . 18°). This will allow the extraction of angular 
distributions starting from Q ss 60 MeV up to and beyond the r\ threshold (see 
Fig. 16). In addition, the fully symmetric WASA detector covers the entire 
azimuthal angular range of 2ir and, thus, is well suited to measure polarization 
observables. The simulations are done for the measurement at Q = 60 MeV 
(studying the development of p- waves) . 

As indicated previously the reaction dd — » an will be identified by detecting 
the a and the decay ir° — ► 77 in coincidence. Here, the high energy loss of 
the a in combination with two neutral hits in the calorimeter provides a very 
efficient first level trigger. For a momenta above 1 GeV/c all particles pass the 
(currently used) Forward Window Counter (FWC) of WASA and are stopped 
in one of the layers of the Forward Trigger Hodoscope (FTH) . Subsequent layers 
(e.g. of the Forward Range Hodoscope, FRH) can be used as veto layers for 



29 



0.6 

O 

^ 0.5 

u 

a. 

0.3 
0.2 
0.1 








0-18° ...*■**" TThrashold 

..-■ : * 




A Resonance A 




Q=60NaV /^V 






0=2° \ 


. i . T . T . . . J i . 1. . i . *. . i . .1 . i . . . i . . . i . . . i . . . i V . . 



o.a 



1.2 



1.4 



1.6 



1.S 



2.2 



2.4 26 2.8 

p. / GeV/c 



Figure 16: Kinematics of the reaction dd — > an for various Q values. The 
transverse momentum of the a versus its longitudinal momentum in laboratory 
is plotted. The lines indicate the azimuthal angle covered by WASA. 
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Figure 17: Energy losses of a particles from the reaction dd — > air . Left: 
Reconstructed angle versus total energy loss. Right: Energy loss in the three 
single layers of the Forward Tracking Hodoscope versus total energy loss. 



particles, which arc not stopped. Angles are reconstructed by means of the 
Forward Proportional Chambers (FPC). 

Fig. 17 (right side) and Fig. 18 show the energy losses in the individual layers 
of the tracking hodoscope for a particles as well as for 3 He and tritons. While 
tritons can already be discriminated by these energy loss patterns, 3 He is quite 
similar to a particles 11 . However, a further reduction of the 3 He content can 
already be done by requiring the correct kinematics as shown in Fig. 17 (left 
side) . 

Together with the measured 7r°, the reaction is kinematically over-constrained. 
Tests on the energy and momentum conservation prevents further particles from 
being missing. While the standard procedure would be an overall kinematic fit, 

11 Protons and deuterons are not shown. Their energy loss is — compared with tritons — 
again smaller and the bands for stopped particles are shifted further to lower energy losses. 



30 



i.OI 0J02 0.03 0.04 0.06 0.0B 0.07 OAS 
Total Energy Loss f GeV 



.01 0.02 0.03 0.04 0.05 OJOB 0.07 0.08 
Total Energy Loss f GaV 



Figure 18: Energy loss in the individual layers versus total energy loss. On 
the left 3 He, on the right tritons. While the energy loss pattern for 3 He is 
rather similar to 4 He, the one for tritons can be used efficiently for particle 
discrimination. 




Figure 19: a kinematics reconstructed from measured ir°. Left: Reconstructed 
polar angle versus measured polar angle. Right: Reconstructed a missing mass. 



Fig. 19 shows — as an example — the correlation between the measured polar 
angle of the a and the one reconstructed from the pion kinematics as well as 
the missing mass with only the pion detected. 

Taking into account the current geometry of WASA, the combined acceptance 
for a and 7r° does not vary much and will be e acc ~ 0.5 for all beam momenta 
from 1.2 GeV/c up to rj threshold. For the following estimates an additional 
factor of 0.5 is considered representing dead-time corrections, detector- and 
analysis efficiencies. 

Assuming a luminosity of L = 10 32 cm~ 2 s _1 = 10~ 4 s _1 pb~ the final count 
rate will be 15 events per week and pb. Since the goal of this experimental 
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program is to extract differential cross sections, the total yield has to be at 
least a few hundred events for each beam energy. Using a total cross section of 
a w 75pb at Q = 60 MeV (i.e. scaling the IUCF result by s-wave phase space) 
leads to beam times of about 1 or 2 weeks per measurement. Consequently, 
although a lower luminosity would still allow to measure one data point within 
a longer (but still reasonable) beam time, it would hardly be possible to carry 
out the full experimental program. This program would cover measurements 
in two different energy ranges, namely at Q « 60 MeV and in the A region 
(Q f» 160 MeV). However, initially one data point at the lowest beam energy 
should be measured. It is aimed to run with polarized beam to disentangle 
the contribution of different partial waves (especially p-waves) as discussed in 
section 2.1.2. Further runs will then depend on the analysis of these data. 
As shown above the reaction dd — > an can be measured with the existing 
WASA setup without any modifications on the detector system itself. The 
electronics has to be upgraded: the trigger logic has to be adopted in order to 
use an efficient energy loss trigger and the data acquisition systems has to be 
replaced to be able to run at high luminosities (see section 4.3.4). 

3.5 «q-/o mixing in pn — >d7r rj and dd — >air°r] 

At a later stage of the experimental program with WASA at COSY the reactions 
pn — ► dir°n and dd — + aTr n will be measured. The main challenges for the 
identification of the isospin violating effects in these reactions will be: 

pn — ► dTT°T]: The measurement of this reaction is similar to the pp experiment 
described in Section 3.2. However, an additional complication comes from 
the fact that deuterium has to be used as an effective neutron target. The 
signal for isospin violation in this reaction would be provided from the 
measurement of the angular forward-backward asymmetry of the (ir°n)i=o 
system. 

dd — ► cxTv°r]: This reaction is only possible due to isospin violation and the 
yield will be rather small. First estimates show that the cross section 
should be of the order of 100 pb [111]. A better prediction can be obtained 
from data on the reaction dd— >aK + K~ which are expected from ANKE 
in winter 2004/05. It is thus suggested to perform these measurements 
(at lower beam energies, see Sect. 3.4) before the dd — > air a rj reaction will 
be studied at maximum COSY momentum (« 3.7 GeV/c). 

3.6 Study of hyperon resonances 

In order to understand the nature of the A(1405) hyperon resonance it is pro- 
posed to study A(1405) production in proton-proton collisions, and in a second 
step also in proton- nucleus collisions. The primary goal of this study is to mea- 
sure the A(1405) production cross section and, in particular, its spectral shape. 
As a by-product also information on the £(1385) resonance is obtained. 
Due to its large acceptance for both charged particles and photons, the WASA 
detector allows to investigate all isospin combinations of E7r and An that are 
populated in hyperon resonance decay subsequent to the production reactions 
pp — > A(1405)if+p and pp — > [E(1385)K] + p. With the use of a deuterium 
target the implementation of a proton spectator detector close to the interaction 
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point would be required, in order to allow the full measurement of all particles 
in the final state (equivalent to pp collisions) also in pn — > A(1405)A'°p and 
pn — > [£(1385)AT]°p reactions. Note that the latter reaction also allows to 
excite the £~(1385) state. 

The photon detection capability of WASA matches particularly well to the study 
of the pp — > A(1405)if + p reaction with the decay channel A(1405) — > £°7r° — > 
(A7X77) which is not populated in S°(1385) decay, and thus should reflect the 
A(1405) spectral distribution in an undisturbed way. Therefore special emphasis 
should be given to this channel, but as discussed in Sect. 2.2.2 also the £ + 7r _ 
and £~7T + channels deserve investigation in order to extract a possible isovector 
contribution. For the observed final state in the decay channel A(1405) — > 
£°7r° — * (A7X77) — > (p7r~7)(77) the expected branching ratio is rather large, 
and amounts to 21%. 

Fig. 20 shows the simulated transverse versus longitudinal momentum distribu- 
tion for the particles in the final state in this reaction induced by 3.6 GeV/c 
protons, assuming a pure phase space distribution of the events, and taking into 
account a A(1405) Lorentz mass distribution of a full width T = 50 MeV/c 2 . 
It demonstrates that both protons in the final state, from the primary vertex 
and from the A decay, have a high probability to be emitted into the acceptance 
of the forward detector, whereas the dominant fraction of the pions will be de- 
tected by the central detector of WASA. The K + mesons populate both regions 
of the phase space seen by the central and by the forward detector. Only a 
small fraction of photons is not covered by the Csl calorimeter acceptance. 
In order to add new information to the experimental knowledge of the spectral 
shapes of the A(1405) and £(1385) resonances, the mass resolution has to be 
clearly better than the widths of these resonances of f=a 50 MeV/c 2 and sw 
40 MeV/c 2 , respectively. The resolution of better than sw 10 MeV/c 2 allowed 
by the WASA detector is considered to be sufficient for this study. Further 
improvement of the resolution is possible if new tracking detectors, such as /it- 
strip silicon arrays, are implemented close to the target. This would allow, in 
addition to the reduction of non-strange background, to make use of additional 
kinematic constraints by a precise measurement of tracks including displaced 
vertex information. 

The optimum proton beam energy is given by the requirement that the kine- 
matic limit be at sufficient distance from the pole of the A(1405) resonance in 
order to guarantee the undisturbed measurement of its spectral distribution. 
At energies too close to threshold it may be difficult to disentangle the A(1405) 
spectral function and the effect of an energy dependent final state interaction. 
A reasonable beam momentum is p v = 3.6 GeV/c, well within the COSY mo- 
mentum range for internal beams, at which the kinematic limit in the A* or £* 
mass is M max = 1.53 GeV, more than 100 MeV/c 2 above the A(1405) pole. 
The expected hyperon resonance production cross sections are of the order 
of 1 fib [191]. With a luminosity L = 10 31 cm _2 s _1 , assuming 10% over- 
all efficiency, one thus expects approximately 18000 counts per day for the 
pp — ► A(1405)-K" + p reaction with a (p7r~7) s o (77)^0 K + p final state discussed 
above. In case the cross section and/or efficiencies should be smaller, the maxi- 
mum design luminosity L = 2-10 32 cm _2 s _1 of WASA gives room for an increase 
of the luminosity. 

The second part of the program related to comparative studies of the in-medium 
properties of the A(1405) and £(1385) hyperon resonances requires the usage 
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Figure 20: Simulated transverse vs. longitudinal momentum distribution for 
the final state particles in the reaction pp — > A(l 405)pK + — ► £ Q 7r°pA" + — ► 
A77r°pA' + — ► p^-K^( b )j( c )^j( d )p( aS > K + (f'. The superscripts denote the corre- 
sponding panels. The ellipses in panels (a) and (f) show the kinematic limits 
with A(1405) production at its pole mass as well as 25 MeV/c 2 below and above. 
The straight lines separate the acceptance of the central and the forward spec- 
trometer of WASA. 



of nuclear targets, and thus an extension of the existing pellet target to an 
operation with gases like nitrogen, argon and xenon which still needs to be 
developed. Furthermore, for these studies additional tracking information close 
to the target allowing hyperon and K s identification will be necessary, since a 
full reconstruction of the final state is not possible with a nuclear target, which 
results in a loss of kinematic constraints. For this purpose the WASA detector 
would have to be extended by the implementation of //-strip silicon detectors in 
close geometry to the interaction point. 

3.7 Rare and very rare r] and rj' decays 

High precision branching ratio for 7/ — > 77. In spite of the large branching 
ratio this experiment will be performed at a later time since it requires technical 
modifications to achieve high precision as discussed below. The radiative decay 
widths of 77(77') — > 77, together with additional information from radiative J/^! 
decays to 7?(t/)7 or from radiative transitions of vector (V) and pseudoscalar (P) 
mesons in V — > Pj and P — ► Vj processes, allow to extract the fundamental 
constants of the pseudoscalar nonet (see e.g. [24]). The unprecedented statistical 
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accuracy, that is expected to be feasible for radiative decays using the WASA 
facility at COSY (see table 1) suggests to increase the precision on branching 
ratios of radiative decays presently available. 

For the decay rj' — > 77, the statistical accuracy of the branching ratio might be 
increased by an order of magnitude in a dedicated run of two to three weeks, 
yielding a statistical error below the percent level. However, to obtain preci- 
sion values for the branching ratio requires to control the systematic error with 
comparable accuracy. A suitable experimental approach is described for the 
decay 77 — > 77 in [192], with a direct measurement of the branching ratio. Thus, 
all systematic uncertainties related to the "normalization", i.e. the simultane- 
ous measurement of a known branching ratio, are eliminated, and, analogously 
to [192] the branching ratio is given by 



analysis and in the electronics to reconstruct and reliably digitize an 77' — > 77 
decay. The precise number of events, background subtraction and (different) 
efficiency corrections only have to be considered for the 77 decay mode. Ex- 
perimentally, the approach requires to trigger on the pp system associated with 
77' production, i.e. on two protons in the forward detector. With an estimate of 
a 50 kHz rate from a trigger on two charged particles in the forward detector 
at luminosities of 10 cm s , and in view of the attainable event rates for 
the planned data acquisition system (see section 4.3.4) in the order of 10 kHz 
a high precision measurement of the 77' — > 77 branching ratio seems feasible 
after all, but will require further improvements of the selectivity of the forward 
detector trigger. Developments have already started for WASA at CELSIUS to 
implement a missing mass trigger for the AE/E technique by implementing en- 
ergy reconstruction and angular correction of the energy loss information on the 
trigger level, and direct branching ratio measurements should become feasible 
after some experience and development with WASA at COSY. 

Double vector meson dominance in 77 — > e + e~e + e~. Experimentally, the 
77 will be tagged by using the missing mass of the two protons detected in the 
forward detector. Tracks of electrons and positrons are measured in the Mini 
Drift Chamber (MDC) inside the super-conducting solenoid, and their energy 
is determined by means of the Csl calorimeter. 

Background from 2 7r° production via 7r°7r° — > e + e~e + e~77 with two undetected 
7s can be effectively removed by kinematic cuts. The 77 decay modes to e + e~7 



(T(ji -» e+e-j)/T tot = 6.0±0.8-10~ 3 [45]) and to 2 7 (T( V -» 77)/^* = 39.43± 



0.26% [45]) contribute via photon conversion in the beam tube. Both channels 
are effectively reduced to a level below 1 % by reconstructing the vertices of 
both e + e~ pairs [193]. 

Based on the QED prediction for the branching ratio of 2.52 • 10 -5 [51] and the 
acceptance of the WASA setup, at a luminosity of 10 32 cm _2 s _1 we expect a 
rate of about 3 • 10 3 events per week. In a 3-months run one will collect a data 
sample of 30000 events which will allow for a high statistics measurement of the 
form factor dependence on the virtual photon masses. In a second phase the 
data for the decay r\ — > e + e~e + e~ can be collected together with the 77 — > e + e~ 
experiments. 




(10) 



where A v i^ 11 denotes the detector acceptance and e v > 



•77 



the efficiencies in the 
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Search for CP violation in the rare decay 77 — ► 7r + 7r e + e . r\ tagging 
as well as electron and positron detection will be performed as described in 
Section 3.1. Electrons are separated from pions using the momentum over en- 
ergy ratio with an accuracy of 5 • 10~ 3 . The fraction of misidentified charged 
leptons has been estimated to be less than 2.5 • 10~ 5 . From the complete four- 
momentum information for all four charged decay products, the angle between 
the 7r + 7r _ and e + e~ production planes can be reconstructed. 
The major source of background arises from the conversion of photons from the 
rj — > 7r + 7r~7 decay in the beam tube. With an overall conversion probability of 
0.28 % the background is important only in the region of small e + e~ invariant 
masses below 40MeV/c 2 . In comparison, background from other sources like 
direct production is negligible. 

The experiment can be carried out together with studies of the rj — ► e + e~e + e~ 
decay channel (see paragraph above) and statistics of nearly 4 • 10 5 r\ — > 
■K + ir~e + e~ events could be obtained, corresponding to an accuracy of 2 • 10~ 3 
for the asymmetry. 

Search for physics beyond the Standard Model in the very rare 
•q — > e + e~ decay. The experiment aims at reaching a sensitivity that allows to 
test whether the 77 decay to the e + e~ mode occurs at a rate compatible with the 
Standard Model prediction of 5 • 10 -9 [43,44], which makes the decay suscep- 
tible to contributions from physics beyond the Standard Model, with a present 
experimental upper limit at 7.7 • 10~ 5 [45]. 

The experimental signature will consist of four charged tracks, with two protons 
in the reaction pp — > ppr/ — > ppe + e~ tagged in the forward detector and the 
e + e~, with an isotropic distribution in the 77 rest frame measured in the central 
detector with magnetic field. In view of the low branching ratio the experiment 
requires high acceptance and high luminosity, charged particle tracking, and 
the possibility to reject events with photons by means of an electromagnetic 
calorimeter, which makes the WASA facility at COSY a unique device for this 
rare decay channel. The most important physics processes that can contribute 
to a possible background are the rj — > je + e~ Dalitz decay, the radiative decay 
mode 77 — > 77, and prompt e + e~ production. The continuous spectrum of e + e~ 
pairs up to the rj mass from the "/e + e~ decay mode (T(ry — > 'ye + e~)/T to t = 
6.0 ± 0.8 ■ 10~ 3 [45]) is effectively reduced to a level of less than 10 % compared 
to the 77 — > e + e~ signal by rejecting photons in the central detector (Fig. 21). 
7 conversion in the material between the interaction point and the plastic 
barrel (PSB) leads to a background contribution from the ?/ — > 77 mode 
(T(r) — > 77)/r tot = 39.43 ± 0.26% [45]), if both 7s convert asymmetrically in 
energy. Most of the relevant material is concentrated by the 1.2 mm Be beam- 
pipe, giving rise to a signal-to-background ratio 12 of 1-2 to 1 in the 77 mass 
range (Fig. 21). A statistically significant signal from the rare 77 — * e + e~ decay 
on the level of the Standard Model prediction, as indicated in figure 21, could 
be obtained in a running time of 2 months. It should be noted, that an ap- 
propriate trigger setting will allow to run the experiment in parallel with other 
77 decay studies, and, that already after one hour of beamtime, the experiment 
can increase the level of sensitivity in comparison to the present upper limit 

12 The background is expected to be reduced further by vertex reconstruction, which is not 
implemented in the simulations. 
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Figure 21: Background contribution from ?y — > 7e + e~ (black), 77 — ► 77 (dark 
grey) and sum of the 77 — > e + e~ signal and the two background sources (light 
grey). The simulation assumes the Standard Model prediction for the branching 
ratio of 5 • 10 -9 for 77 — > e + e~ and the PDG values for the background channels, 
scaled to a total of 10 10 77 events produced. 



by an order of magnitude. Furthermore, it should be emphasized, that these 
investigations are only feasible due to the optimization of the WASA setup with 
respect to the amount of material between the interaction point and the cen- 
tral detector (see section 4.2), to fully exploit the advantages of a windowless 
internal target. 

However, the theoretical and experimental uncertainties for the background 
from prompt e + e~ production in the 77 mass range are large. Calculations [194] 
based on the approach described in [195] lead to a signal-to-background ratio for 
the decay 7/ — > e + e~ of 1 to 30 with respect to prompt production. Thus, in order 
to have a reliable estimate for the prompt production background, it is suggested 
to prepare the 77 — > e + e~ search by measuring the prompt contribution pp — > 
ppe + e~ first. 

At the same time, since WASA allows for exclusive measurements, this prepara- 
tory study can shed light on both the relative strength of different sources of 
e + e~ pairs, and meson and baryon isobar form factors. The first aspect is inti- 
mately motivated by the expectation that dileptons carry information from the 
early stage of heavy ion collisions. Discrepancies between theoretical calcula- 
tions (see e.g. [196,197]) and the presently available data for dilepton production 
in proton-proton and proton-deuteron collisions measured with rather small ge- 
ometrical acceptance at the DLS spectrometer [198] — especially in the mass 
range between 0.2 and 0.5GeV/c 2 — could be verified experimentally using 
WASA at COSY. Assuming a luminosity of 10 32 cm~ 2 s _1 reasonable statistics 
of about 10 4 events can be achieved in five days of data taking, i.e. a program 
of four kinetic beam energies ranging from 1.04 to 2.5 GeV to overlap both with 
the DLS data and theoretical calculations is feasible within three weeks, and 
will provide the necessary information for a reliable background estimate from 
prompt production concerning the search for the rare 77 — > e + e~ decay. 



37 



C violation in 77 — ► 7r°e + e~. The experiment aims at a sensitivity down to 
the Standard Model prediction for the branching ratio of 0.2— 1.3 • 10~ 8 [52-54], 
to search for new C-violating effects in the electromagnetic interaction. With 
a value of T(rj — *■ 7r°e + e~)/r tot < 4 • 10~ 5 [45] the present upper limit is three 
orders of magnitude weaker than the theoretical estimate. 

Monte Carlo simulations have been carried out to investigate the feasibility of 
the experiment and background conditions. In conclusion, a search sensitivity 
of 10~ 9 is possible considering the full background. A sample of ten events from 
C-conserving mechanisms for the rare decay can be expected after a period 
of three weeks to five months, depending on the actual value of the branching 
ratio. 

A C-violating admixture at the level of 10 -3 would already double the 77 — > 
7r°e + e _ decay rate. However, at present this effect is small compared to the un- 
certainty in theoretical estimates for the C conserving decay mechanism, which is 
due to insufficient knowledge of value and structure of the 77 — ► 7r°77 amplitude. 
Consequently, an improved precision in the 77 — » ir°jj decay rate is important 
for the interpretation of the experimental result of the proposed study of the 
ir°e + e~ decay mode, and a high statistics study of 77 — > 7r°77 should be carried 
out simultaneously. 



Measurement of the if transition form factor. With 77' tagging by means 
of two detected protons in the forward detector and the reconstruction of the 
decay system of the Dalitz conversion decay rj' — > 77* — > 7e + e~ in the cen- 
tral detector, the WASA setup is well suited for measuring the 77' transition 
form factor. With the tracking information from the Mini Drift Chamber, the 
momenta of electron and positron are determined. Thus, the differential cross 
section 
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[F(q 2 )r (11) 

pointlike 



with respect to the four-momentum transfer squared q 2 of the time-like virtual 
photon, which is equal to the invariant mass squared of the e + e~ pair, can be 
measured, and the transition form factor F(q 2 ) can be extracted. The form 
factor is usually fitted with the simplest pole-type formula 

F(q 2 ) = — ^ « 1 + £ (12) 

with the slope parameter A. 

At present, only an upper limit exists for the r )e + e~ decay mode of the 77' 
(T(?/ -> "fe+e-)/T tot < 9- 1CT 4 ). With an estimated branching ratio of 3- 1(T 4 , 
as discussed in [55] , we expect 45 fully reconstructed if Dalitz conversion decays 
per day. Following the discussion for the rf — ► ~/Tr + ir~ decay in section 3.1 we 
only consider the radiative decay if — > 77 with subsequent conversion of one 
7 as a possibly important source of background. With fits to the 77' signal in 
each bin of the e + e~ invariant mass scale, the differential cross section can be 
extracted background free, unless there is a background process for the eta' — ► 
7e + e~ transition. However, since for 77 decays the 77 mode was found to be 
negligible above very low invariant lepton pair masses [37] , and since the ratios 
T(j](j]') — > 7e + e _ )/r(77(77') — > 77) arc approximately the same for 77 and 77', and 
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since the photon conversion probability is not changing drastically, it can safely 
be neglected as a background contribution. 

The beam time estimate depends on the desired accuracy for the slope parameter 
A in eq. 12. From the analogous discussion of the r\ transition form factor 
in [37] , a running time of two months would allow a determination of the slope 
parameter with a 10 % error. A further improvement by a factor of two in 
accuracy would require a nine months period of data taking. It should be noted, 
that due to the trigger condition of two charged and one neutral particle in the 
central detector data taking could be done in parallel with other rf decay studies. 
For example, the specific trigger requirement for the rf Dalitz conversion decay 
exactly matches the trigger setup for the 7r + 7r~7 mode discussed in section 3.1 
in terms of particle multiplicities. 
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4 Experimental facility 



4.1 Cooler Synchrotron COSY 

COSY is a cooler synchrotron and storage ring operated at the IKP of FZJ. The 
accelerator complex (see Fig. 22) comprises an isochronous cyclotron (JULIC), 
used as an injector, a race track shaped cooler synchrotron with a circumference 
of 184 m, and internal and external target stations [199]. COSY delivers beams 
of polarized and unpolarized protons and deuterons in the momentum range 
between 0.3 GeV/c and 3.7 GeV/c. The ring can be filled with up to 10 11 par- 
ticles leading to typical luminosities of 10 31 cm _2 s~ 1 when using an internal 
cluster target. Beams can be phase-space cooled by means of electron cooling 
at injection energy as well as stochastic cooling at high energies. Typical beam 
preparation times, including injection, accumulation and acceleration, are of the 
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Figure 22: Floorplan of the accelerator complex. The WASA detector will be 
installed at the place of the old cavity in front of the electron cooler. 
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order of a few seconds, while the beam lifetime with a cluster target is between 
several minutes and an hour. 

Currently, four internal experiments (ANKE, COSY-11, EDDA, PISA) and 
three external detector systems (BIG KARL, JESSICA, TOF) are operated by 
large international collaborations. On the average COSY is running for more 
than 7000 hours per year. Typically, it delivers beams for experiments with a 
reliability of 94%. 

Implementation of WASA at COSY 

The Cooler Synchrotron has two 40 m long straight sections joining the arcs. 
After removal of the old RF cavity at the beginning of 2004 space for the imple- 
mentation of WASA at an internal target position in front of the electron cooler 
is available (see Fig. 22). This location has a number of advantages: most im- 
portantly it minimizes both the interference with the existing COSY detectors 
and the modifications needed to the WASA detector and the accelerator. 
Operating WASA with the pellet target with an effective thickness of 2.5 • 
10 15 cm -2 average luminosities of 10 32 cm _2 s _1 and beam lifetimes of a couple 
of minutes are expected [200] . Compared to the present operation at CELSIUS 
the experimental conditions will be improved due to the different accelerator 
characteristics: fast magnet ramping, dispersion-free target position, stochastic 
cooling and a smooth microscopic time structure of the beam. 

4.2 The WASA detector 

The 47r detector facility WASA [201,202] was designed for studies of production 
and decays of light mesons at CELSIUS. Pions and eta mesons are produced 
in proton-proton and proton-dcuteron interactions. The highest beam-proton 
kinetic energy reachable at CELSIUS is 1.5 GeV. 

The pellet-target system is integrated in the setup and it provides small spheres 
of frozen hydrogen and deuterium as internal targets. This allows high luminos- 
ity and high detection coverage for meson decay products like photons, electrons 
and charged pions. 

WASA consists of a forward part for measurements of charged target-recoil 
particles and scattered projectiles, and a central part designed for measurements 
of the meson decay products. The forward part consists of eleven planes of 
plastic scintillators and of proportional counter drift tubes. The central part 
consists of an electromagnetic calorimeter of CsI(Na) crystals surrounding a 
superconducting solenoid. Inside of the solenoid a cylindrical chamber of drift 
tubes and a barrel of plastic scintillators are placed. A vertical cross section of 
the WASA detector is shown in Fig. 23. 

4.2.1 Pellet target 

The pellet target system was a unique development for the CELSIUS/WASA 
experiment [203,204]. The main components of the system are shown in Fig. 24. 
The heart of the setup is the pellet generator where a jet of liquid hydrogen is 
broken up into droplets with a diameter around 35 //m by a vibrating nozzle. 
The droplets freeze by evaporation in a droplet chamber and form a beam of 
pellets that enter a 7 cm long vacuum-injection capillary. After collimation, the 
pellets are directed through a thin 2 m long pipe into the scattering chamber 
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Figure 23: Cross section of the WASA detector. The central detector built 
around the interaction point (at the left) is surrounded by an iron yoke. The 
layers of the forward detector are visible on the right-hand side. The individual 
components are described in the text. 



pellet 
beam lube 

CELSIUS 
beam lube 



pump station L 




pellet generator 
differential pumping 

scattering chamber 



n n 



pump station 
cryogen icbeamdump" ~~~~ - 

Figure 24: Layout of the pellet target system. 



and further down to a pellet beam dump. The inner diameter of the pipe is 5 
mm at its entrance to the scattering chamber. This arrangement provides the 
necessary space to put the Ait detection system around the interaction region. 
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Pellet target thicknesses of up to 3 • 10 15 atoms/cm 2 give acceptable half- lives 
of the circulating ion beam as well as acceptable vacuum conditions. Some of 
the parameters of the pellet target at the present stage of operation are listed 
in table 3. The pellet target system operates regularly with pellets of normal 
hydrogen and deuterium. 



Pellet target parameters 


Present performance 


Pellet diameter 
Pellet frequency 
Pellet-pellet distance 
Effective target thickness 
Beam diameter 


25 - 35 fim 
5-12 kHz 
9-20 mm 
> 10 15 atoms/cm 2 
2-4 mm 



Table 3: Present performance of the pellet target system. 



4.2.2 Forward detector 

The forward detector (FD) is designed mainly for detection and identification 
of scattered projectiles and charged recoil particles like protons, deuterons and 
He nuclei in n and r\ production reactions. Also neutrons and charged pions can 
be measured. All FD plastic scintillators may supply information for the first 
level trigger logic. This part of the setup has already been used in a previous 
experiment at CELSIUS, WASA/PROMICE, which is described more in detail 
in [205]. A summary of the most important features of the forward detector 
(FD) is given in table 4. The individual components are described in some 
detail in the following. 



Forward detector 


Total number of scintillator elements 


280 


Scattering angle coverage 


3° - 17° 


Scattering angle resolution 


0.2° 


Amount of sensitive material [g/ cm 2 ] 


50 


[radiation lengths] 


PS 1 


[nuclear interaction lengths] 


ps 0.6 


Thickness of vacuum window (st. steel) [mm] 


PS 0.4 


Maximum kinetic energy (T s top) for stopping: 




7r ± /proton/dcuteron/alpha [McV] 


170/300/400/900 


Time resolution 


< 3 ns 


Energy resolution for: 




stopped particles 


ps 3% 


particles with T stop < T < 2T st0 p 


4-8% 


Particle identification 


AE-E 



Table 4: Some features of the Forward Detector. 



The Forward Window Counters (FWC) 

The FWC is the first detector layer in the FD (along the beam direction) and 
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consists of 12, 5 mm thick plastic scintillators (Fig. 25). It is mounted tightly 
on the paraboloidal stainless steel vacuum window. Therefore, the elements are 
inclined by approximately 10° with respect to the plane perpendicular to the 
beam direction. The FWC signals are used in the first level trigger logic to 
reduce the background caused by particles scattered in the downstream beam 
pipe and in the flange at the entrance to the FD. The signals arc also used to 
select He ejectiles on the trigger level. 




Figure 25: Schematic view of the FWC. 



The Forward Proportional Chambers (FPC) 

Immediately downstream of the FWC, there is a tracking device. It is composed 
of 4 modules, each with 4 staggered layers of 122 proportional drift tubes (so 
called straws) of 8 mm diameter (Fig. 26). The modules are rotated by 45° 
with respect to each other (in the plane perpendicular to the beam axis). They 
are used for accurate reconstruction of track coordinates and provide precise 
angular information of the particles originating from the target region. 

The Forward Trigger Hodoscope (FTH) 

The FTH, consisting of three layers of 5mm thick plastic scintillators, is placed 
next to the FPC. There are 24 Archimedian spiral shaped elements in the first 
two planes and 48 radial elements in the third. 

The FTH is mainly used in the first level trigger logic. However, the special 
geometry, combining all three layers, results in a pixel structure, which is useful 
for resolving multi-hit ambiguities. In Fig. 27 the structure of the FTH is shown 
with hits from two passing charged particles. 
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view of one module of the FPC. 




Figure 27: Schematic view of the three forward trigger hodoscope layers. 
The Forward Range Hodoscope (FRH) 

Behind the FTH, the four layers of the FRH are positioned (Fig. 28). Each plane 
is made of 24, 11 cm thick, plastic scintillator modules. The FRH, together 
with FTH, is used for energy determination of charged particles and for particle 
identification by AE-E technique. Fig. 29 shows how protons and deuterons can 
be identified. 

The identity and initial kinetic energy of a charged particle is reconstructed 
from the pattern of deposited energy in the different detector planes. Fig. 30 
shows how the initial kinetic energy of a proton is related to the deposited en- 



45 



Figure 28: Schematic view of the four forward range hodoscopc layers. Dimen- 
sions are given in mm. 
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Figure 29: Example of particle identification in the forward detector for data 
collected at 640 MeV. Energy deposited in the I s * layer of FRH is plotted versus 
total energy deposited in FTH and FRH. The lines indicate cuts to be used for 
selecting protons (lower) and deuterons (upper). 



ergy in the forward detector scintillators. Even for particles that are stopped 
in a detector, the total deposited energy is different from the initial kinetic en- 
ergy. This is because some of the energy is lost in inactive material between 
the detectors elements, in the scattering chamber windows, etc. Fig. 30 also 
shows that the variation of the deposited energy is strong enough to be use- 
ful for energy reconstruction also for high energy particles not stopping in the 
detector material. For protons, this can be used in the kinetic energy range 
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Figure 30: Energy absorbed in the forward detector scintillator material for 
protons as a function of their initial kinetic energy. The upper curve gives the 
total energy deposited in the plastic scintillator planes and the lower curves 
show the energy distributed among the different FRH planes. 



300 McV to 800 MeV and for dcutcrons in a similar energy range starting from 
400 McV. Identification of punch-through particles can be done using either the 
veto hodoscope (FVH) or the AE information from the last FRH planes. 

The Forward Range Interleaving Hodoscope (FRI) 

Between the third and fourth layers of the FRH there are two interleaving 
layers of 5.2 mm thick plastic scintillator bars (Fig. 31). Each layer has 32 bars, 
oriented horizontally in one and vertically in the other. 

The main purpose of this addition to the FRH is to provide a two-dimensional 
position sensitivity inside the FRH necessary for measurement of scattering 
angles for neutrons. The efficiency for detection of neutrons of a few hundred 
MeV kinetic energy by the FRH is around 35 %. In addition, the FRI can help 
in vertex reconstruction and to discriminate against background tracks due to 
secondary interactions in the beam pipe and other structural material. 
The FRI was recently commissioned; more information about the design and 
performance of the FRI can be found in [206] . 

The Forward Veto Hodoscope (FVH) 

The last detector layer of the FD is a wall of plastic scintillators (Fig. 32). It 
consists of 12 horizontal plastic scintillator bars, equipped with photomultiplicrs 
on both sides. The hit position along a bar may be reconstructed from signal 
time information. In the first level trigger the signals are used for rejection (or 
selection) of particles punching through the FRH. 

The Forward Absorber (FRA) 

Optionally, a passive absorber layer made of iron can be introduced between 
the last layer of the FRH and the FVH. The thickness of the absorber can be 
chosen from 5 mm up to 100 mm. The absorber has been used for stopping the 
protons from the pp — > ppr\ reaction at a beam proton energy around 1360 MeV. 
In this case the faster protons from elastic scattering and from pion production 
penetrate the FRA and induce signals in the FVH which can be used for veto 
purposes in the trigger. 
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Figure 31: Schematic view of the FRI detector (upper picture) and its two 
planes with orthogonally oriented scintillator bars (lower picture). 




Figure 32: Schematic view of the Forward Veto Hodoscope. 
4.2.3 Central detector 

The central detector (CD) is built around the interaction point and is de- 
signed mainly for detection and identification of the decay products of tt° and 
•q mesons: photons, electrons and charged pions. It consists of an inner drift 
chamber (MDC), a solenoid (SCS) providing magnetic field for momentum mea- 
surements, thin plastic scintillators in a cylinder geometry (PSB) and a Csl 
calorimeter (SEC). The amount of structural material is kept minimal to re- 
duce the disturbances on the particles. The beam pipe is made of 1.2 mm thick 
beryllium and the total thickness of the solenoid corresponds to 0.18 radiation 
lengths only. 

For the design the main requirements were the following: 

• it has to handle high particle fluxes at luminosities around 10 32 cm~ 2 s -1 . 



48 



• it should be able to measure photons with energies from a few MeV up to 
800 MeV. 

• it should be able to measure, in a magnetic field of about 1 T, the mo- 
menta of electrons and positrons in the range p ss 20 — 600 MeV/c with 
an accuracy cr p /p < 2%. 

The momenta of heavier charged particles can also be measured in a similar 

momentum ranges, but with lower accuracy: 

for pions and muons with p rts 100 — 600 MeV/c, er p /p < 4% 

for protons with p w 200 - 800 MeV/c, <r p /p < 6% 

The main components of the Central Detector, shown in Fig. 23, are presented 
below in some detail. 

The Superconducting Solenoid - (SCS) 

The SCS provides an axial magnetic field for the momentum measurements for 
the tracks measured by the MDC. It also protects the CD from low-energy delta 
electrons produced in the interactions of beam particles with the pellets. The 
wall thickness of the SCS is minimized in order to allow high accuracy of the 
energy measurements in the calorimeter. The return path for the magnetic flux 
is provided by a yoke made out of 5 tons of soft iron with low carbon content. 
The yoke shields the readout electronics from the magnetic field and serves also 
as support for the calorimeter crystals. The main parameters of the SCS are 
given in table 5. 

In order to map the magnetic field inside the volume enclosed by the SCS, the 
magnetic field strength inside the MDC was measured with Hall probes and, in 
addition, the field distribution was calculated with simulation programs. The 
calculated values were fitted to the measured ones with an error of ±1% of 
Btotai [207]. The magnetic flux density distribution inside of the iron yoke 
calculated for a current of 667 A is given in Fig. 33. The SCS is described in 
detail in the Ph.D. thesis of Roger Ruber [208]. 



Superconducting coil 


Inner/outer radius [mm] 


267.8 / 288.8 


Superconductor (stabilizer) 


NbTi/Cu (pure Al) 


Total winding length 


465 mm 


Maximum central magnetic flux density, B c 


1.3 T 


Field uniformity in the MDC 


1.22 T ±20% 


Cooling 


Liquid He, 4.5°K 


Cryostat 


Material 


Aluminium 


Inner / outer radius [mm] 


245 / 325 


Overall length [mm] 


555 


SCS wall thickness (coil+cryostat) [radl] 


0.18 



Table 5: Main parameters of the superconducting coil and its cryostat. 
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Figure 33: Calculated distribution of the magnetic flux density for a coil 
current of 667 A [207]. Contour maxima are indicated by lines marked 
A-H, where: A=0.10T, B=0.25T, C=0.50T, D=0.75T, E=1.00T, F=1.20T, 
G=1.30T, H=1.50T. 

The Mini Drift Chamber - (MDC) 

The MDC is placed around the beam pipe and is used for determination of 
particle momenta and reaction vertex. It is a cylindrical chamber covering 
scattering angles from 24° to 159°. For large angle scattered protons from elastic 
proton-proton scattering, a vertex resolution (er) of 0.2 mm perpendicular and 
3 mm along the beam axis can be reached. A detailed description of the MDC 
can be found in the Ph.D. thesis of Marek Jacewicz [209]. 

The MDC consists of 1738 drift tubes, so called straws, arranged in 17 cylindrical 
layers. The diameter of the straws in the 5 inner most layers is 4 mm, 6 mm in 
the next 6 layers and 8 mm in the outer 6 layers. The straws are made of a thin 
(25 (im) mylar foil coated with 0.1 fim aluminum on the inner side only. In the 
center of each straw there is a 20 fim diameter sensing wire made of gold plated 
tungsten (W(Re)), stretched with a tension of 40 g. The wires are aligned with 
a precision of ±20 /xm. 

This design was chosen in order to cope with the expected high particle flux 
allowing a maximum deposited energy of approximately 70 ^m^s ^ or ^ ne s t raw s 
exposed most at the inner part of the chamber. 

The layers are located between radii of 41 and 203 mm. The straws in nine 
layers are parallel to the beam axis (z-axis) and the other eight layers have 
small skew angles (6°-9°) with respect to the z-axis. These "stereo" layers form 
a hyperboloidal shape. 

The straws in the five inner layers are divided unequally by the center of the 
pellet pipe, while the other layers are symmetrical. The straws in each layer are 
inter-spaced by small gaps in order to prevent the mechanical deformation by 
neighboring tubes. 

The MDC is fitted inside a cylindrical cover made of 1 mm Al-Be and is placed 
inside the solenoid (Fig. 34). 

The straws in each (half) layer are mounted between «5 mm thick Al-Be end- 
plates. The layers are assembled around the Be beam-pipe (Fig. 35) and the 
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Figure 34: (Left) The fully assembled MDC inside the Al-Be cylinder. (Right) 
The MDC surrounded by PSB elements. 




Figure 35: (Left) Drift tubes secured in the end-plates. Note the "stereo" layers 
interleaved with parallel layers. (Right) Be beam-pipe. 

attached pipe for the pellets. The beam-pipe has a diameter of 60 mm and a 
wall thickness of 1.2 mm. The design drawing of the Be beam pipe is shown in 
Fig. 36. 
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Figure 36: Beryllium beam-pipe including the crossing pellet target tube. 



The Plastic Scintillator Barrel - (PSB) 

The PSB is located inside the SCS coil and surrounds the MDC. It provides fast 
signals for the first level trigger logic and, together with the mini drift chamber 
and the Csl calorimeter, it is employed for charged particle identification by the 
AE-p and AE-E methods and serves as a veto for 7 identification. 

The performance of the PSB has been studied using proton-proton elastic scat- 
tering events. Fig. 37 (left plot) shows the result of a Monte Carlo simulation 
of the angular dependence on the energy deposited in the PSB. Fig. 37 (right 
plots) shows typical experimental spectra after a correction for nonuniform sig- 
nal response has been applied. 

In the initial experiments the momentum and energy resolution allowed reason- 
able discrimination between pions and protons, which is illustrated in Fig. 38. 
For high energy charged particles also SEC information is available and can be 
used for the identification. 

The PSB consists of a cylindrical part and two end caps and contains in total 
146 elements shaped as strips of 8mm thickness. In the cylindrical part there 
are 48(+2) elements of 550 mm length and 38 mm width, forming 2 layers with 
a small (on average 6 mm) overlap between neighboring elements to avoid that 
particles pass without registration. The end caps with an outer diameter of 
approximately 42 cm in the backward and 51 cm in the forward part contain 48 
"cake-piece" shaped elements each. The front end cap is flat while the rear cap 
forms a conical surface. Both end caps have a central hole for the beam pipe 
(19 cm diameter in the forward and 12 cm diameter in the backward part). The 
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Figure 37: (Left) Angular dependence of the deposited energy in the PSB for 
simulated elastically scattered protons. The energy deposition increases with 
increasing polar angle (corresponding to a decrease of the kinetic energy of 
the proton), until particles begin to stop in the plastic scintillator material (at 
around 9=77°). (Right) Experimental spectra corrected for light attenuation 
for four of the PSB central elements. 




Figure 38: Example of particle identification in the central detector for a raw 
data sample collected at 1360 MeV. Energy deposited in the PSB is given as 
a function of signed momentum from the MDC. The regions for protons and 
pions are marked. 

PSB as modeled in the detector simulation program is shown in Fig. 39. One 
sector of the PSB is shown in Fig. 40 and 41. 

Each scintillator is glued to an acrylic light guide coupled to the photomultiplicr 
tube (PMT). The PMTs are placed outside of the iron yoke to shield them from 
the magnetic field. For this purpose, approximately 50 cm long light guides are 
used. 

The Scintillator Electromagnetic Calorimeter - (SEC) 

The CD calorimeter SEC is able to measure photons, electrons and positrons 
with energies up to 800 MeV. The energy threshold for detection of photons 
is about 2 MeV. SEC consists of 1012 sodium-doped Csl scintillating crystals 
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Figure 39: Forward, central and backward parts of the PSB. In the central part, 
the gaps for the pellet target pipe are visible. 

Type "C": Forward Type "D": Backward 




Figure 40: (Left) Layout of one section of the PSB detector. A denotes the 
rectangular counters of the barrel wall and B and C are trapezoidal elements 
in the forward and in the backward caps respectively. D are bended light 
guides. (Middle) Two shapes of the trapezoidal forward elements with dimen- 
sions marked in mm. (Right) Shape and dimensions in mm of the trapezoidal 
backward element. 

placed between the super-conducting solenoid and the iron yoke. The scattering 
angles covered by the SEC arc between 20° and 169°. 

The crystals are shaped as truncated pyramids and are placed in 24 layers along 
the beam (Fig. 23 and 43). The lengths of the crystals vary from 30 cm (16.2 
radiation lengths) in the central part to 25 cm in the forward and 20 cm in the 
backward part. Fig. 42 shows the angular coverage together with the thickness 
of SEC. As a measure of the anticipated photon fluxes, the center of mass (CM) 
system solid angle vs. the laboratory (LAB) scattering angle is shown for some 
experimental conditions at WASA. 

A planar map of the crystals is given in Fig. 44. One can clearly distinguish the 
three different main parts of the calorimeter: the forward, central and backward 
parts. The forward part consists of 4 layers with 36 elements each. It covers 
scattering angles from nearly 20° to 36°. The central part consists of 17 layers 
each having 48 elements, and covers scattering angles from 36° to 150°. The 
backward part consists of three layers. Two layers have 24 elements and one 
layer closest to the CELSIUS beam pipe has only 12 elements. The small spaces 
between the forward-central and central-backward parts are occupied by PSB 
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Figure 41: Cross section of the PSB scintillators with dimensions marked in 
mm. No structural material is shown. 
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Figure 42: The angular coverage of the SEC. The CM system solid angle vs. the 
LAB scattering angle is shown for pp and pd interactions at 1500 McV and 
895 McV. 



light guides and mechanical support for the solenoid (back end only). The 
calorimeter covers nearly 360° in <i> angle. Holes for the pellet pipe (2+2 crystals) 
and the solenoid chimney (4 crystals) are not shown in the figure. Some design 
parameters of the calorimeter are given in table 6. 

The SEC is composed of sodium-doped Csl scintillating crystals. This type of 
scintillator material provides a large light yield, has short radiation length and 
good mechanical properties. CsI(Na) was chosen instead of the more commonly 
used CsI(Tl) scintillators for the following reasons [210,211]: 
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Figure 43: Schematic view of the SEC. It consists of the forward part (shadowed 
area to the left), the central part (not shadowed area in the middle) and the 
backward part (shadowed area to the right). The beam is coming from the right 
side. 
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Figure 44: Planar map of the SEC. 



• Its emission peak at 420 nm matches well the bi- alkali Sll photocatode of 
the selected PM tubes, giving good photon statistics and sufficiently fast 
response. 
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Scintillator Electromagnetic Calorimeter 


Amount of sensitive material 


135 g/cm 2 


[radiation lengths] 


« 16 


[nuclear interaction length] 


ps 0.8 


Geometric acceptance: 


96% 


polar angle 


ps 20° - 169° 


azimuth angle 


ps 0° - 360° 


Max kinetic energy for stopping 




71"^ /proton/deuteron 


190/400/500 


Scattering angle resolution 


« 5°(FWHM) 


Time resolution 




charged particles 


5 ns(FWHM) 


photons 


«40 ns(FWHM) 


Energy resolution 




charged particles 


ps 3%(FWHM) 


photons 


ps 8%(FWHM) 



Table 6: SEC design parameters. 



• Its shorter scintillation decay time is preferable in high-rate applications. 

• CsI(Na) gives much less afterglow than CsI(Tl). 

• CsI(Na) seems more resistant against radiation damage. When irradiated 
by a proton beam corresponding to 10 years of operation a test crystal 
did not show any visible change in its structure. The CsI(Tl) test crystal, 
on the contrary, lost its transparency. 

The crystals are connected by plastic light guides, 120 mm to 180 mm long, 
with the photomultiplicrs placed the outside of the iron yoke. In Fig. 45, a fully 
equipped single calorimeter module consisting of a Csl crystal, a light guide, a 
PM tube and a high voltage unit, enclosed inside a special housing, is shown. 
The SEC and its performance is described in more detail in the Ph.D. thesis of 
Inken Koch [183]. 

The light pulser system (LPS) 

The LPS delivers reference light pulses via light fibers to all scintillation counters 
in order to monitor their gain during the experiment. Since both organic and 
inorganic scintillators are used, two types of light sources were designed. A 
xenon flash tube from Hamamatsu is used for the Csl elements of the calorimeter 
and three LED-based light sources for all plastic scintillators. From those four 
sources the light signals are transmitted to individual elements via a network of 
light fibers [212]. 
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Figure 45: A fully equipped module with Csl crystal, light guide, PM tube and 
housing. 

4.3 Modifications 

4.3.1 Forward detector (FD) 

The FD was designed mainly for measurements of protons of kinetic energies up 
to 500 MeV. An energy resolution cte/E around 3% is obtained at the highest 
energy. The track coordinates are measured by the FPC straw chambers. At 
CELSIUS, three modules with 4 layers of straws each, are installed. To improve 
track coordinate measurements, energy measurements and particle identification 
for the higher energies anticipated at COSY, the following modifications are 
planned: 

• The first FRH layer will be removed. It gives the few centimeters of extra 
space, needed for the installation of one more (existing) FPC module. By 
moving the remaining FRH layers a few centimeters upstream, one gets 
a better energy measurement performance at the largest scattering angles 
(17°) covered by the FD. 

• Repair of FTH. The FTH and FRH have served at CELSIUS for twelve 
years and signs of severe aging effects have appeared. Most detectors still 
work well, but in one of the spiral planes of FTH some elements do not 
give acceptable signals and the whole plane should be replaced. 

• Two new FRH planes should be installed to improve particle energy re- 
construction. This would mean that protons of about 350 MeV kinetic 
energy, instead of 300 MeV at present, would be stopped in the FD (fig- 
ure 46). The energy resolution at high energy would consequently also be 
improved. A cte/E around 10% could be reached for protons at 900 MeV. 

• Downstream of the FRH, we consider the installation of a 10 cm thick 
water Cherenkov detector, that is composed of w2 m long and 10 cm 
wide bars. For protons it would have a range of sensitivity from 500 
MeV (/3=0.75) to 1200 MeV (/3=0.9), which well covers the region where 
the energy determination based on the FRH depositions deteriorates. The 
combination of the FRH and the Cherenkov information will improve both 
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the particle identification power and the kinetic energy resolution. As an 
example, the Cherenkov light output is plotted as a function of kinetic 
energy for protons and pions from the reaction pp — > pprj' where rf — > 
rjir + 7r~ at a beam momentum of 3.3 GeV/c (figure 47). The proton kinetic 
energy energy resolution is about 5% for protons above 500 MeV. Studies 
are underway to obtain more quantitative information on the performance 
of the proposed Ccrcnkov detector. 
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Figure 46: Energy absorbed in the (upgraded) FRH planes for protons (upper 
curve) and pions (lower curve) as a function of their kinetic energy at the target. 



4.3.2 Central detector (CD) 

The CD was originally designed for measurements of gammas with energy up 
to 600 MeV and of electrons and positrons with momenta up to 400 MeV/c. 
The charged particle momenta are determined by the track curvature in a mag- 
netic field of about 1 Tesla. The maximum lever arm is 160 mm for these 
measurements and the design accuracy in the position measurements in each of 
the 17 MDC layers is 100 /itm. The accuracy of the field map should be better 
than 1%. The design momentum resolution for electrons and positrons, scat- 
tered at large angles is < 2% (er p /p). With the present implementation of the 
MDC in the CELSIUS/WASA experiment an optimal resolution, a p /p «3% for 
e , w4% for tt^~ and «8% for protons could be obtained. In real experiments, 
the momentum resolution has been evaluated for and protons and cr p /p 



59 




EKIN [MEV] 

Figure 47: Chcrcnkov light output for charged particles from the reaction pp — > 
pp?7' at 3.3 GeV/c as a function of kinetic energy. 

values of «4.5% and wl2% were obtained [209]. To improve the momentum 
measurements for work at higher energies the following actions arc planned: 

• An increase of the maximum magnetic field to 1.3 Tesla. This puts great 
demands on some parts of the magnet system, since there will be a 70% 
increase in the stored energy and in the mechanical forces involved, but it 
would not require any major hardware upgrades [208]. However, it needs 
careful tests and "training" of the magnet by expert personnel. 

• The MDC electronics have to be upgraded in order to improve the effi- 
ciency for detection of pions and electrons and to reach the design time 
resolution of w0.5 ns. The extent of this task depends on the actual elec- 
trical noise environment at the WASA site at COSY. 

The design energy resolution for photons measured in the SE calorimeter is 
about 5 % (<je/E). At present an effective energy resolution of 7-10 % is 
obtained [183], and the main limitation is a non-linear signal response. This 
is the product of some different effects and the situation will be improved by 
modification of PM bases, by usage of new QDCs and, if neccessary, by changing 
some 150 PM tubes. 

One source of background in CELSIUS/ WAS A is interactions in the restgas in 
the scattering chamber. At present there are three big cryopumps at WASA, 
two at the forward cone of the scattering chamber and one at the backward 
cone. The vacuum in the scattering chamber will be improved by adding one 
more big cryopump in the backward direction. This needs a careful redesign of 
the backward part of WASA. 

Another source of background is due to scattering events originating from the 
vacuum chamber walls. Most of these events are caused by beam halo particles. 
It is suggested to put a ring-shaped detector inside of the beampipe just at the 
entrance to WASA. Signals from this detector should be used in the trigger logic 
to "veto" events caused by beam halo. Further prestudies will be done to design 
such a detector and evaluate its anticipated performance in more detail. 
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4.3.3 Cost estimates 



Cost estimate 


kEuro 


FRH upgrade with two additional layers 


200 


FTH repair / replacement of one layer 


50 


FD Cerenkov detector 


100 


SEC new PM bases / change of PM tubes 


50 


MDC system refurbishments 


130 


Additional upstream vacuum pump 


40 


Upstream internal veto detector 


30 



Table 7: Cost estimates for some planned detector modifications. 



4.3.4 DAQ- and Trigger-System 

The WASA DAQ-System was designed more than 10 years ago and implemented 
with electronics from the same period. Most of the digitizing modules were 
delivered by the company LeCroy, which is no longer active on the market of 
nuclear electronics. Furthermore, many parts of the detector electronics have a 
very low level of integration, leading to a high cabling effort at the detector side, 
especially for the straws in the MDC and the Forward Detector. In addition, 
spare parts of the used components cannot be acquired anymore. 
Due to the progress in technology the readout speed of the existing system is 
much slower than possible using a modern, state-of-the-art system. Since the 
event rate at COSY must be significantly higher than the one currently reached 
at CELSIUS, this becomes a serious problem. Considering also the maintenance 
problems and the age of the digitizing modules, it is obvious that successful 
operation of WASA at COSY can only be achieved by a major upgrade of the 
existing DAQ system. 

Currently, all experiments running at COSY are using the same DAQ system 
following a uniform approach developed, implemented and maintained by per- 
manent staff of FZ Julich. In view of the time scale for "WASA at COSY" the 
upgrade of the WASA DAQ should be based on the third generation of the DAQ 
at COSY. The main advantages are: 

• it implements the desired functionality, including a modern synchroniza- 
tion technique, avoiding the - inherently slow - sequential readout of the 
digitizing modules, 

• it is almost completely developed, thus fitting into the limited time scale, 

• it has a good on-site support by the original developers, which is essential 
for a smooth and successful operation, 

• and it is compatible to all existing DAQ systems at COSY. 
The current DAQ system 

The existing DAQ system of WASA at CELSIUS is mainly based on four FAST- 
BUS crates containing the digitizing modules (see Fig. 48). The readout is done 
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by PCs connected to the FASTBUS crates via a proprietary parallel link based 
on RS485. Since there is no synchronization system, all digitizing modules have 
to be read out sequentially when a trigger occurs. The readout of the (about) 
1500 PMTs is done with LeCroy LRS1881 QDCs. The time information from the 
straws (2000 channels from the FPC and 1700 channels from the MDC) and the 
w 500 fast plastic scintillators is digitized using LeCroy LRS1876 TDCs [213]. 
The trigger system of WASA at CELSIUS is based on dedicated hardware that 
has been developed by TSL. Due to the high background the trigger system 
is still evolving to reach higher performance. The first level trigger uses the 
discriminator outputs of the fast plastic scintillators. It is based on multiplicity 
and generates the gates for the QDCs and the stop for the TDCs. The second 
level trigger uses the discriminator outputs from the calorimeter. It is based on 
cluster multiplicity and energy deposition and generates the FastClear for the 
digitizing modules. 




Figure 48: Structure of the existing WASA DAQ and trigger system [213,214]. 
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The new DAQ system 

Starting with the first experiments at COSY permanent staff from ZEL was 
responsible for a common DAQ system and its further development. In the 
meantime two evolutionary and elementary steps have been made to adapt it to 
the rapidly changing technologies [215]. Now, the implementation of a third gen- 
eration of DAQ has been started aiming for the highest possible event rates [216]. 
Therefore, state-of-the-art FPGA technologies are combined with fast commu- 
nication paths to achieve system latencies as low as possible. Commonly used in 
all three generations is a DAQ software called EMS [215], which guarantees the 
software compatibility of all implemented systems. EMS follows a well proven 
client /server architecture and has been developed during a period of more than 
10 years (corresponding to a manpower effort of nearly 30 man- years). 
As in other DAQ systems (e.g. at CERN) this third generation of DAQ at COSY 
is based on proprietary solutions in order to guarantee implementations opti- 
mized for speed. However, also CAMAC and VME can be integrated despite of 
the resulting loss of performance. The whole system is designed to run exclu- 
sively in a "common stop mode" to avoid many expensive and hard to handle 
delay units. Here, "common stop mode" means that the trigger is available 
some 100 ns to 1 /is after the signal has been digitized. Therefore, the acquisi- 
tion boards must run in a self-triggering mode digitizing all interesting signals 
and storing them together with a time stamp. When the - also time-stamped 
- trigger arrives an FPGA selects the digitized signals inside a predefined time 
interval. 

An essential key component required by this operation mode is a Synchroniza- 
tion System (SYS) to reliably control and synchronize the fast event flow [217]. 
The SYS delivers a global time base relative to the trigger time. The control of 
the synchronization system as well as the readout of the digitization modules is 
done by a farm of embedded industrial PCs. On this farm, EMS-servers con- 
trolling the front end modules are implemented. As a future extension also the 
cancellation of events (software trigger) is considered. 

The third generation of DAQ at COSY is based on a bus controller (refer to 
Fig. 49) connected by a fiber optical link to the PC farm and the synchronization 
system. The acquisition modules are connected by a high speed LVDS bus. This 
bus can be implemented as a backplane or as a cable bus covering a distance of 
15 meters and serving for control as well as for transmission of event data. The 
physical layer of the bus complies to the SCSI bus with a maximum transmission 
speed of 80 or 160 Mbytes/s [218]. In general, acquisition, trigger processing, 
data reduction, buffer management, and sub-event building can be done by 
means of FPGAs — either on the acquisition boards or in the bus controller. 
All these functions are realized without the intervention of a general purpose 
processor and, thus, lead to the envisaged performance. 

As shown in Fig. 50 a new computing layer and a new layer of digitization 
electronics will be implemented. The readout will be done by a farm of 15 PCs 
connected via an optical link (a common development of ZEL and the company 
SIS, with a physical layer identical to Gigabit Ethernet). The overall experiment 
control and storage is connected via Gigabit Ethernet. The whole system runs 
under control of the EMS software as used at all other experiments at COSY. 
The digitization layer consists of crates with readout controllers, QDC and TDC 
modules and the synchronization system. The event buffering capabilities of the 
digitization modules in combination with the synchronization system allow the 
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Figure 49: Structure of the third generation DAQ system at COSY. 



operation in common stop mode, thus avoiding delay lines. It is expected that 
with the new system the event rate will be increased at least by a factor of 
10. Because of performance requirements and costs the use of commercially 
available digitizing modules is not possible. 

The existing 64 channel TDC module (being used for the proportional chambers 
of ANKE and the straws at TOF) can be used without modifications for the 
straws at WASA, provided the discriminator electronics will be changed to a 
CMP16 based system. It is a 6U board for the LVDS bus based on the Fl 
chip system developed by ZEL compatible with the third generation of DAQ at 
COSY [219]. 

Optionally, also a VME board developed at TSL may be used, although some 
modifications and a further test phase are necessary [220]. This board has 64 
channels, too, and is based on the CERN TDC32. The decision which TDC to 
use will have only minor impact on the overall cost estimate and project timing 
and will be done within the next year. 

For the calorimeter (pulse duration about 2 fis) the QDC will be implemented by 
sampling the analog signal with a 100 MHz FlashADC and subsequent integra- 
tion in an FPGA. Considering connector requirements and layout restrictions, 
it is planned to have 16 channels on one 6U base board conforming to the 
standards of the LVDS bus system developed by ZEL. 
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Figure 50: Structure of the new DAQ and trigger system of WASA at COSY. 



The latter technique is not possible for the signals from the fast plastic scin- 
tillators, because the pulse duration is only « 20 ns. Here, additional shaping 
or even the use of a QAC chip is required. The technical issues are intensively 
discussed with people from the University of Giessen (TAPS collaboration), the 
KFKI Budapest and the TSL. Independent of the final decision, a mezzanine 
board will be manufactured fitting to the FlashADC board described above and 
performing the analog preprocessing for 16 channels. 

The time stamps required by the free-running mode will be provided by the 
FPGA on the motherboard for the "slow" QDC using the global time informa- 
tion and by additional discriminators and TDCs for the "fast" QDC. 
As indicated in Fig. 50, the discriminators and the trigger electronics from the 
existing system shall remain. Only with regard to the straws it is intended to 
replace the discrete preamplifier and discriminator electronics with new modules 
based on the CMP16 chip. The main advantage will be the increased sensitivity 
and the compact layout, which allows much shorter cable lengths between straws 
and preamplifier, thus reducing noise. 

A further improvement of the trigger system can be achieved by implementing 
a third level trigger in the context of the new DAQ system. While this is not 
possible in the limited time scale of about one year, it should be considered as 
a future extension. At the moment two basic hardware concepts are discussed: 

• Implementation as a processor- or FPGA-farm that operates directly on 
the DAQ data stream. 
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Electronics for MDC and FPC based on the CMP16 ASIC 






3700 channels 




Digitizing modules 






TDCs (MDC, FPC) 


3700 channels 




QDCs (slow) 


1000 channels 




QDCs (fast) 


400 channels 




add. TDCs 


400 channels 




System components 






17 Crates (incl. crate controllers, interfaces) 




Industrial PC farm (embedded) 


30 k€ 


Communication equipment 


20 k€ 


Synchronization system 


70 k€ 


Miscellaneous 






Racks, infrastructure test and measuring equipment 


80 k€ 


Experiment control, 


visualization and storage 


20 k€ 


Additional cables and connectors 


70 k€ 


External developments 


50 k€ 


6% risk reserve 




60 k€ 


Total 




1000 k€ 



Table 8: Cost estimate for the development of a new data acquisition system 
for WASA at COSY. 



• Splitting the data from the digitizing modules into two independent 
branches, one for the trigger and one going to the DAQ system [221]. 

In a further future-oriented development project it is planned to integrate the 
embedded PC functionality in the LVDS bus controller itself using "system on 
chip" (SOC) technology with a general purpose FPGA and processor on the 
same chip. However, it is quite uncertain, whether this project will be finished 
next year. Therefore it has not been considered for the DAQ system of WASA 
at COSY. 

Cost estimate 

The cost estimate shown in table 8 does not include manpower, since all de- 
velopments concerning the new DAQ system and improvements of the existing 
WASA electronics can be made by existing permanent staff. One exception 
could be the development of the analog part of the fast QDCs. Most of the ex- 
isting cables should be reused. However, additional shielded cables are required 
in some cases and all connectors will be replaced. 
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5 Project plan 



General remarks 

The project "WASA at COSY" deals with the relocation of an existing detector 
system — the Wide Angle Shower Array (WASA) from The Svedberg Lab- 
oratory (TSL) at Uppsala university (Sweden) to the Research Center Jiilich 
(Forschungszentrum Jiilich, FZJ, Germany), more specifically to an internal 
target position of COSY (Cooler Synchrotron) at IKP. WASA at COSY will be 
the experimental activity with the highest priority at the IKP in Jiilich for the 
next several years. 

The approach of adapting an existing detector has a number of advantages 
compared to building a new one: It is faster, provided the detector is available 
on short notice. Given the necessary lead-time for preparations at COSY, the 
date that WASA will be available (end of 2005) is regarded as optimal. It is 
cost-effective if the detector is "in good shape" which is the case for WASA. 
However, the movement and setting up of a detector as complex as WASA is not 
free of charge — even neglecting basic costs such as dismounting, transport, and 
installation. We anticipate the following investments for necessary repairs and 
changes to the detector: Adaptation to the higher COSY energy, in particular 
the forward detectors of WASA. Exchange of electronics for which no spare 
parts exist and no replacement can be purchased any longer. Replacement of 
the data acquisition system to improve data taking capabilities and to adapt 
the WASA-DAQ to COSY standards. These will allow one to carry out the 
initial experiments described above. 

Manpower 

Well over one hundred scientists have committed themselves to actively persue 
the physics program laid out in this proposal. This includes members of the 
current WASA at CELSIUS collaboration as well as a wide representation from 
the COSY users. As the host institute, IKP will have a special role in providing 
the scientific and technical infrastructure for WASA at COSY. Consequently, 
IKP will redirect significant manpower (about 10 F.T.E.) to this project, which 
implies a reduction of the other activities at IKP. 

Finances 

The financing of this project is not yet finalized. In view of the high priority 
WASA at COSY has for IKP about one half of the needed investments summa- 
rized below could be contributed from the running budget of IKP. This would 
require a one-time reduction of the number of COSY running hours by 50% for 
one year in either 2005 or 2006. The remaining investment must be provided 
by other sources in either Sweden, the EU, the German BMBF, or the Research 
Center Jiilich. 



Dismounting, transfer, set-up 100 k€ 

Preparations at COSY 300 k€ 

Changes at WASA ("adaptation") 600 k€ 

New trigger- and readout electronics 1000 k€ 

Sum: 2000 k€ 
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Time schedule 



June 2005 Termination of CELSIUS operations at TSL 

Oct. 2005 Finish dismounting of WASA at CELSIUS 

Nov. 2005 Transfer of WASA to Jiilich 

June 2006 Finish set-up of WASA at COSY 

July 2006 Start of commissioning 

Jan. 2007 Start of experiments 



(>8 
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